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Chapitre 1 : Vocabulaire thermodynamique 


Introduction 


D'origine grec, le mot thermodynamique est composé de deux parties : «thermo» qui veut 


dire chaleur et «dynamique» qui signifie travail ou mouvement. Le mot composé veut dire 


mouvement produit à partir de la chaleur. La thermodynamique est apparue en tant que 


science a partir du XVIII siècle avec la création de la machine à vapeur. Son but est d’étudier 


les différentes formes d’énergies et les possibilités de conversion entre elles : énergie 


thermique <> énergie mécanique, chimique.....etc. 


La thermodynamique compte plusieurs branches, à savoir : 


Thermodynamique classique (considère les systèmes sous leurs aspects macroscopiques 
en utilisant des grandeurs mesurables telles que la pression P, la température T et le 
volume V) 

Thermodynamique chimique (étude des chaleurs des réactions et des propriétés 
thermodynamiques des équilibres chimiques) 

Thermodynamique statistique (considère les systèmes sous leurs aspects 
microscopiques et s’appuie sur la mécanique statistique) 

Thermodynamique technique (application de la thermodynamique a létude des 


machines thermiques, frigorifiques,..etc) 


Quelques dates historiques 


1747 : le physiologiste A. HALLER a attribué la constance de la température, du corps 
humain au frottement accompagnant la circulation du sang. 

1765 : J. WATT a créé la premiere machine a vapeur pratiquement sans fondement 
théorique. 

1824 : S. CARNOT a découvert le deuxième principe de la thermodynamique, qui a été 
complété plus tard (en 1848) par W. THOMSON et R. CLAUSIUS (en 1850). 

1842 : R. MAYER a découvert le premier de la thermodynamique, qui a été confirmé par 
J. JOULE en 1845. 

1869 : M. MASSIEU, J.W. GIBBS (1875), H.L.F. HELMOTZ (1882), et P. DUHEM 
(1886) ont découvert les fonctions caractéristiques et les potentielles thermodynamiques, 
et ont initié l’application de la thermodynamique à la chimie. 


1877 : L. BOLTZMANN a introduit la thermodynamique statistique. 
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e 1906: W.NERST a présenté le troisième principe de la thermodynamique. 
e 1931 : L. ONSAGER a introduit la thermodynamique des processus irréversibles. 


1.2 Système thermodynamique 
Un système est une quantité de matière (solide, liquide ou gaz) séparée du milieu extérieur 
par une surface appelée frontière à travers laquelle se fait un échange d’énergie (chaleur Q, 


travail W) et de matière. 


muheu extérieur 





SYSTEME <-> énergie - matière 


suface 


Figure 1.1 : Système thermodynamique 


La nature de la surface (ou frontière) du système thermodynamique sert de base à la 
classification des différents systèmes : 
e Systeme ouvert (fig.1.2a): échange d’énergie (sous forme de travail W ou de chaleur Q) 
et de matière avec le milieu extérieur. 
e Systeme fermé (fig.1.2b): échange d’énergie seulement avec le milieu l’extérieur. 
e Systeme isolé (fig.1.2c) : aucun échange d’énergie ni de matière avec le milieu extérieur. 
e Systeme adiabatique (fig.1d): aucun échange de chaleur avec l’extérieur (surface 
imperméable à la chaleur) 


e Surface diathermane : surface (ou paroi) permettant le passage de la chaleur. 


F== Hy pip n Ca a kae p a eN, 

















compresseur 
d'air 





(a) Sys. ouvert (b) Sys. ferme (c) Sys. isolé (d) Sys. adiabatique 


Figure 1.2 : Schéma montrant les différents types de systèmes thermodynamiques 
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1.3 Etat d’un système 


L’état d’un système est defini par les valeurs d’un certain nombre de grandeurs physiques 
mesurables dites variables (ou paramètres) d’état. Exp: volume, pression, température, 


viscosité,...etc. 


1.3.1 Variables d’état 

Les variables qui définissent l’état d’un système peuvent être classées en deux catégories : 

e Variables intensives : ce sont des grandeurs indépendantes de la quantité de matière 
(ne sont pas additives) et définies ponctuellement, exp: la pression (P, [kPa]), la 
temperature (7, [K]), la masse volumique (p, [kg/m° J), potentiel chimique (u, 
[kJ/kmol]),..etc. 

e Variables extensives : ce sont des grandeurs proportionnelles à la quantité de matière 
(sont additives), exp : la masse (m, [kg]), le volume (V, [m°]), le nombre de mole (n, 
[kmol]),..etc. 

Par exemple, supposons deux systèmes identiques placés ensemble. Il est clair que les 
variables extensives (m, V et n) vont doubler alors que les variables intensives (P, T et p) ne 


vont pas changer. 


1.3.2 Fonctions d’état 

Ce sont des relations mathématiques qui lient les variables d’état (intensives ou extensives) 
d’un système thermodynamique, exp: l’énergie interne (U, [kJ]), Venthalpie (H, [kJ]), 
l’entropie (S, [kJ/K]),..etc. Les fonctions d’état sont extensives et leurs variation, au cours 
d’une transformation, dépend uniquement de l’état initial et final du système et non du 


chemin suivi (la manière dont le système passe de l’état initial a l’état final). 


Remarque : une fonction d'état (extensive) rapportée à l'unité de matière (masse, volume, .. 
etc) devient intensive. Par exemple, l'énergie interne (U, [kJ]) est extensive mais l'énergie 


interne massique (u=U/m, [kJ/kg|) est intensive. 


1.4 Equilibre thermodynamique 

Un système est en équilibre thermodynamique lorsque toutes ses variables d’état intensives 
demeurent uniformes (identiques en tout point du système) et constantes dans le temps. 
L’équilibre thermodynamique résulte de la conjonction de trois équilibres : thermique, 
mécanique et chimique. 

e Equilibre thermique (égalité des températures) : la température, en tout point du 


système, uniforme et constante. 
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Equilibre thermique 





Etat fmal 





Etat mitial 


Figure 1.3 : Equilibre thermique d’un système 


e Equilibre mécanique (égalité des pressions): la pression, en tout point du système, 


uniforme et constante. 


Equilibre mécanique 





Etat initial Etat final 


Figure 1.4 : Equilibre mécanique d’un système 


e Equilibre chimique (égalité du potentiel chimique) : la composition chimique du 


système uniforme et constante. 


PIE 


E | fi Te 





Intermédiaire 
Figure 1.5 : Equilibre chimique après diffusion de l’encre dans l’eau 
1.5 Transformation thermodynamique 


C’est le passage d’un état d'équilibre initial (1) à un état d’équilibre final (2) au cours 


duquel au moins une variable d’état est modifiée. 
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Figure 1.6 : Représentation graphique 
d’une transformation thermodynamique 





1.5.1 Transformation irréversible 

Si au cours d’une transformation, qui fait passer le système d’un état d’équilibre (1) à un 
état d’équilibre (2), les états successifs intermédiaires ne sont pas des états d'équilibre, la 
transformation est dite irréversible. Notons que toutes les transformations réelles sont 


irréversibles. 


Transformation rapide 





équilibre 
final 


Figure 1.7 : Compression irréversible par dépôt brusque de la masse M 


Les causes d’irréversibilité sont : 
- les frottements de tous ordres (mécanique, visqueux..) ; 
- le transfert de chaleur (corps placé dans un four...) ; 
- les phénomènes de diffusion (encre dans l’eau, fuite d’un gaz comprimé vers 
l’extérieur..) ; 


- les réactions chimiques. 


1.5.2 Transformation quasi-statique 
C’est une transformation qui fait passer un système d’un état d’équilibre à un autre par une 
succession d’états d’équilibres. De telles transformations ne peuvent avoir lieu que de 


manière infiniment lente (puisque l’état d’équilibre doit être réalisé a tout moment). 
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1.5.3 Transformation réversible 

Une transformation réversible est une transformation quasi-statique qui peut évoluer en 
sens inverse (transformation inversable) et passe exactement par les mémes états d’équilibre 
intermédiaires aussi bien pour le système que pour le milieu extérieur. Dans la pratique, une 
transformation peut approcher la réversibilité si elle est effectuée aussi lentement que possible 
et si les frottements sont négligeables, et si les écarts de température (comme tous les 


gradients des variables intensives) restent faibles 





Transformation trés lente 


Figure 1.8 : Compression réversible par dépôt progressif de la masse M 


1.5.3 Transformations particulières 

e Transformation isotherme : transformation à température constante. 

e Transformation isochore : transformation à volume constant. 

e Transformation isobare : transformation a pression constante. 

e Transformation adiabatique : transformation qui s’effectue sans échange de chaleur 
avec le milieu extérieur. 

e Transformation fermée (cyclique) : transformation à la fin de laquelle le système 
revient a son état initial. 


Ces transformations peuvent être réalisées de façon réversible comme de façon irréversible. 


1.6 Pression et température 

1.6.1 Pression et échelle de pression 
C’est la mesure de la force exercée par unité de surface d’un système (P=F/S). Elle est 
due aux collisions (chocs) des molécules avec les surfaces du système. En 


thermodynamique, on utilise différents types de pressions : 


e Pression atmosphérique (Pam): pression du milieu ambiant mesurée à l’aide d’un 
baromètre (Fig1.9) 
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AP = Patm — Ps = pg 





Figure 1.9 : Principe du baromètre a mercure 


e Pression absolue (P.p;): toute pression mesurée a partir du 0 bar (le vide) 
e Pression effective (Pep): c’est l’écart entre la pression absolue et atmosphérique, 


mesuré à l’aide d’un manomètre. 


Peff = Pabs — Patm 


Faim 





j Y pE 
. liquide 


Figure 1.10 : Manométres pour mesure des pressions 


e Pression du vide (Pride): c’est l’écart entre la pression atmosphérique et la pression 


absolue, mesuré à l’aide d’un indicateur de vide. 


Pode mn ~ Favs 


. L'unité SI de la pression est le pascal, abréviation (Pa) : 


1 Pa=1N/m?=1kg/ms? , 1 Pa=10° bar , latm= 101325 Pa 
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P ihs = P i 


Paff 
P P p 


atin ann atin 


Pi = P aa 


vide 





absolu 


Figure 1.11 : Schéma montrant la relation entre les différentes pressions 


1.6.2 Température et échelle de température 
Du point de vue microscopique, la température est une mesure de l’agitation moléculaire 

(énergie cinétique moyenne) de la substance : plus le mouvement des molécules est important 

plus la température est élevée. Du point de vue macroscopique, la température est une valeur 
numérique qui caractérise le degré d’échauffement de la substance et permet de prédire le 
sens du transfert de chaleur. 

e Echelle Celsius (°C): échelle utilisée dans le monde entier, elle fixe la température de 
congélation de l’eau à 0°C et sa température d’ébullition à 100°C. 

e Echelle Fahrenheit (°F) : utilisée dans les pays anglo-saxons, en particulier les USA. La 
température de congélation de l’eau est fixée à 32°F et sa température d’ébullition à 
2127P. 

e Echelle Kelvin (K) : c’est une échelle de température thermodynamique (indépendante 
des propriètés de la matière) qui repose sur un seul point de référence, 1l s’agit du point 
triple! de l’eau (Tip = 273,16 K) 

e Echelle Rankine (R) : échelle absolue utilisée dans les pays anglo-saxons. 

T(K) = TCC) + 273,15 
T(R) =T(°F)459,67T(R) = 1,8 x T(K)T(°F) = 1,8 x TCC) + 32 
AT(K) =AT(°C) et AT(R) = AT(°F) 


t Le point triple est un état où les trois phases de la matière (solide, liquide et gaz) coexistent 
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K 
i 
vaporisation — + 


de l'eau + 











373.15 


_ 


point tnple 


de l'eau 1l | 
fusion cli 


la glace L 
Figure 1.12 : Comparaison entre l’échelle 


Kelvin et Celsius 


Kelvin 


=- 


Absolute zero —— 


1.7 Travail et chaleur 
Le travail et la chaleur sont des formes d’énergies” en transite et non contenues par le 
système. On les appelle aussi fonction de parcours car ils dépendent du chemin suivi (et non 


de l’état initial et final). L’unité de ces énergies est le joule (J) 


1.7.1 Travail (W) 

Dans ce cours de thermodynamique, on s'intéresse essentiellement au travail mécanique 
des forces extérieures (dit travail des forces de pression ou travail en vase clos : piston...) et 
au travail mécanique technique (travail fait ou reçu par une machine, on l’appelle aussi travail 
utile)’. 

e Travail des forces de pression 
Le travail élémentaire óW des forces extérieures Fu le long d’un parcours dl est : 
OW = Eu x dl 


En intégrant entre l’état initial (1) et (2), on obtient le travail total : 


(2) 
W = | Foxe X di 
(1) 
Appliquons cette relation au cas de la compression d’un gaz 
(2) (2) Və 
W = f Fo X dl = f Poce XSXdl=- | Pere x dV 
(1) (1) Vi 


* L'énergie cinétique, l’énergie potentielle et l’énergie interne sont d’autres formes d’énergies. Ce sont des 
propriétés du système 

? Il y a d’autres formes de travaux, tel que : le travail électrique (dû aux forces électriques), travail dans le champ 
de gravitation (dû à la force de gravitation), travail élastique (forces élastiques), travail magnétique ...etc 
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dV, désigne la variation élémentaire du volume du système. Le signe (-) dans cette 
relation vient du fait que d/ et dV varient dans des sens contraires. W est donc positif lors 
d’une compression (énergie reçue par le gaz ) et négatif lors d’une détente (énergie 
fournie par le gaz). Cette relation est générale, elle est valable quelque soit la 


transformation (compression, détente, réversible, irréversible, ...). 





Figure 1.13 : Travail des forces de pression lors de la compression d’un gaz 


v àP constante: W = —P,,,(V> — V1) 


F | : n < . yV.: 

v à P variable et transformation irréversible (Pext #Poaz) : Wire = — i Pt X AV 
\ i | y. 

v à P variable et transformation réversible (Pext=Peaz) : Wréey = — l Pyaz X dV 


Sur le diagramme de Clapeyron (P, V), le travail (W) mis en jeu lors d’une transformation 
réversibl (1)-(2) correspond à laire entre la courbe qui désigne cette transformation et l’axe 
des abscisses (V). On remarque sur la figure ci-dessus que le travail est une fonction de 


parcours puisque à deux chemins différents deux surfaces différentes. 


transformation 





— 1 du 
Figurel.14 : Travail des forces de pression (W) sur le diagramme de Clapeyron 
v Transformation isochore (V=cst) : W= 0 


v Transformation cyclique : le travail correspond à l’aire du cycle dans le diagramme de 


Clapeyron (Fig.1.15)0W 
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Figure 1.15 : Travail d’une transformation cyclique 


e Travail technique 
Le travail techique est donne par la relation suivante : 
(2) 
W, = f V x dP 
(1) 
Sur le diagramme de Clapeyron (P, V), le travail (W,) mis en jeu lors d’une transformation 
réversible (1)-(2) correspond à l’aire entre la courbe qui désigne cette transformation et l’axe 


des abscisses (P) 


P 






transformation 
réversible 


dP 


l 





Figure 1.16 : Travail technique (Wp sur le diagramme de Clapeyron 


1.7.2 Chaleur (Q) 


La chaleur est une énergie qui apparait lorsque deux corps de température différente sont 
mis en contact. Elle passe du corps chaud au corps froid, JAMAIS l’inverse. On rappelle 
toujours qu'un corps NE POSSEDE PAS une quantité de chaleur déterminée, mais il 
l'échange en fonction des corps avec lesquels il rentre en contact et en fonction du type de 
transformation. 


Dans le cas d’un gaz" qui subit une transformation réversible allant de (1) a (2): 


* Pour les liquides et les solides on écrit simplement 6Q = mc dT (m la masse, c la chaleur spécifique et dT 
l'écart de température). 
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(2) (2) 
Q = f 6Q = f T dS 


(1) (1) 
S, étant l’entropie du système. 


e Lorsque l’échange de chaleur a pour effet de faire varier la température d’un système, on 
parle de chaleur sensible 

e Lorsque la chaleur échangée ne provoque pas de variation de température mais un 
changement de phase (exp : passage de la glace (état solide) à l’eau (état liquide)), on 


parle de chaleur latente (chaleur latente de fusion, chaleur latente de vaporisation) 
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Chapitre 2 : Le premier principe de la thermodynamique 


2.1 Convention de signe pour un système thermodynamique 

On rappelle que la convention sur le signe de l’énergie, en thermodynamique, est telle que 
l’énergie reçue par le système est positive (> 0), et négative (<Q) si elle est cédée par le 
système 





2.2 Les différentes énergies d’un système 


Energie totale 
(Er) 












Energie macroscopique 


y Ÿ 


Energie cinétique 


Energie microscopique 
(ou interne U) 


Energie cinétique Energie potentielle 


Energie potentielle 
(Ec) 


(Ep) 





barycentrique propre d’ interaction 
(Ec) (E p) 





Figure 2.1 : Différentes énergies d’un système 


a) Energie macroscopique (ou externe) 


C’est l'énergie que possède l’ensemble du système par rapport a un référentiel externe 
exp : énergie potentielle et cinétique. 


> l’ énergie potentielle (Ep) est liée à l’ élévation du système dans le champ gravitationnel 


du systè | ae 
masse ike] ss accélération élévation du centre 
terrestre [m/s?] de gravité [m] 
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> l énergie cinétique (Eç) est liée au déplacement du système par rapport à un référentiel fixe 


Ee=m C*/2 [kJ] 


vitesse de déplacement 
du systeme [m/s] 


b) Energie microscopique (ou interne) 


C’est l’énergie liée a la structure moléculaire du système et au degré d’activité des 
molécules (mouvements de translation, rotation, vibration et interaction entre molécules). 
Cette énergie est notée U (énergie interne qui est une fonction d'état! extensive) et elle ne 
dépend d’aucun référentiel externe du système. 


U = Ep + Ec 


énergie potentielle énergie cinétique 
d'interaction barycentrique 


c) Energie totale 


L’énergie totale (Er) d’un système est donc la somme de son énergie macroscopique 
(externe) et microscopique (interne) : 


Er = Ep + Eç +U 
La variation de l’ énergie totale d’un système est de ce fait égale a: 


AE, = AEp + AE; + AU 
ou sous forme différentielle : 


dEr = dEp + dEe + dU 
Remarque 


L’expression de l’énergie totale donnée ci-dessus concerne un système thermomécanique 
(considéré dans ce cours). Pour d’autres types de système, il peut y avoir d’autres formes 
d’énergie (propre au système), telles que: l’énergie électrique (cas de condensateur), 
nucléaire, magnétique, ...etc. 


2.3 Le premier principe de la thermodynamique pour un système fermé 


Lors d’une transformation (réversible ou non) entre deux états d'équilibre I et 2, toute 
variation de l'énergie totale d’un système résulte d’un échange d'énergie sous quelque forme 
que ce soit avec le milieu extérieur. 


' La variation de U dépend uniquement de l’état initial et final de la transformation et non du chemin suivi 
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Q12 + W> 69 t = ÂËr., = AEP,» + AE¢,5 T AU,> (2.1) 
énergies échangées 


variation de l'énergie 
avec l'extérieur 


totale du système 


avec 

O2 énergie thermique 

W2 énergie mécanique (travail de forces de pression) 
C12 énergie électrique 


AEp,, Variation de l'énergie potentielle macroscopique (ou externe) 
AEc,, Variation de l'énergie cinétique macroscopique (ou externe) 


AU;> variation de l’énergie interne 


Comme dans beaucoup de cas, nous avons affaire a des systèmes fermés au repos ou 
stationnaires (AËÇ,, = 0) et dont l’énergie potentielle extérieure est souvent constantes 


(AEp,, = 0), la relation (2.1) se simplifie. Elle devient : 
Qi2 + Wiz + e12 ++ = AET,, = AU 42 


Si de plus nous considérons que les seules formes d’énergies échangées avec le milieu 
extérieur sont de nature mécanique (W) et/ou thermique (Q), cette expression se réduit à : 


Qi2 + Wiz = AE7,, = AU; = U2 — Uy 


Nous avons la une deuxieme forme (utilisée le plus souvent) du premier principe de la 
thermodynamique que nous énonçons de la foçon suivante : 


La somme (algébrique) des énergies thermique et mécanique échangées entre un systéme et le 
milieu extérieur, au cours d’une transformation quelconque (réversible ou non), est égale à la 
variation de son énergie interne 


W= 8k 
W= Sk] 


LAU= (85418) 


= Ski 





O=-12K 


Figure 2.2 : Illustration du premier principe de la thermodynamique 
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Si on évolue entre les états d'équilibre (1) et (2) de deux façons différentes, l’une révérsible 
(a) et l’autre irréversible (b) (voir Fig. 2.3), on peut tout à fait écrire : 


(Q12)a + (Wi2)a = (Q12)p + (Wi2)b = AUs2 = U2 — Us 
S S 

énergies échangées énergies échangées 

de façon réversible de façon irréversible 


(Q12)a # (Qi2)p 
-== Transformation (Wa + (W)p 


~ irréversible 


` b) 


i 









Transformation 
réversible 


(a) 


Figure 2.3 : Représentation d’une transformation réversible (a) et irréversible (b) 


° 212 = 7 . 2 
Pour une transformation élémentaire on écrit : 


6Q + 6W = dE, = dU 


2.3.1 Application du premier principe à quelques transformations particulières 
> Transformation cyclique 


Supposons que dans un premier temps on fournit un certain travail W, au système ci- 
dessous (Fig. 2.4), en faisant tourner un agitateur. Dans un deuxieme temps, on laisse le 
systeme revenir a son état initial par échange d’une quantité de chaleur Q avec le milieu 
extérieur. L’ensemble des deux transformations constitue donc un cycle. 





(a) Q 


(b) 
| 


Figure 2.4 : Application du premier principe de la thermodynamique à un cycle 


L’ application du premier principe de la thermodynamique, nous permet d’écrire : 


2 nites Pe i a oe ae aad : 
On utilise "6" pour désigner la variation infinitésimale (ou élémentaire) d’une grandeur de parcours et "d" pour 
la variation infinitésimale d’une grandeur d’état 
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Q +W =AU = Urinal — Üinitiar = Q (2.2) 


(Étant donné que l’état initial est confondu à l’état final et que U est une fonction d’état = 


Urinal = Vinitiat) 





Vv 


Figure 2.5 : Représentation d’une transformation cyclique 


La relation (2.2) traduit le principe d'équivalence”, qui dit que : 


Si, au cours d’une transformation cyclique, un système quelconque ne peut échanger avec le 
milieu extérieur que du travail et de la chaleur, la somme du travail et de la chaleur échangés 


par le système est nulle 


> Transformation adiabatique (Q = 0) 


Dans ce cas le travail mécanique (W) échangé correspond a la variation de l’énergie 


interne (Fig. 2.6a) : 
W = AU ou óW = dU 


> Transformation isochore (V=Cste) 


Pour un système isochore (parois fixes), la chaleur (Q) échangée avec le milieu extérieur 
correspond à la variation de |’ énergie interne (Fig. 2.6b) : 


Q = AU ou 6Q = dU 


> Transformation isobare (P=Cste) 


Lorsque la transformation est isobare et réversible (Fig.2.6c), la chaleur (Q) échangée avec le 


milieu extérieur correspond à la variation de l’enthalpie : 


Q+W=Q-PAV=AU = 


3 ; ; | es ; Ad > E aah a ee 

Historiquement, le premier principe de la thermodynamique a été appelé principe d’équivalence (ou de 
conservation d’énergie) car il stipulait l’équivalence entre deux formes d’énergie: W (mécanique) et Q 
(thermique) 
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Q = AU + PAV = AH 
Pour une transformation infinitésimale : 


5Q = dH 


> Systeme isolé (Q = 0 et W = 0) 
Pour un système isolé (Fig. 2.6d), la variation de l’énergie interne est nulle : 


0 = AU ou O=dU > U=Cste 


(a) Transformation (b) Transformation | (c) Transformation (d) Système isolé 
adiabatique isochore isobare 


Parois fixes 


Figure 2.6 : Illustration du premier principe de la thermodynamique 
pour des systèmes particuliers 





2.4 Le premier principe de la thermodynamique pour un système ouvert 


Lorsque le système est ouvert, il y a lieu de tenir compte de l’énergie accompagnant la 
matière qui traverse le système (Fig. 2.7) par les différentes canalisations d’entrées (i) et de 
sorties (7). Cette énergie se présente (dans les études habituelles de thermoénergétique) sous 


Su-Yu 


les formes suivantes : 


- énergie interne 


entrée sortie 
2 
Se CF Ci 
- energie cinetique = Mi > 7 My 
entrée sortie 
- énergie potentielle = Mi 9Z; — Mj Zj 
entrée sortie 
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- travail de transvasement” = > PV; — > P;V; 


entrée sortie 
m; etm; désignent, respectivement, la masse de fluide correspondant au flux de matière à 


l'entrée et sortie du système; U; , Ci , Zi , Pi, Vi (et U;,C;,z;,P;,V;) sont l’énergie interne, la 


vitesse, l’altitude, la pression et le volume du fluide lorsque celui-ci traverse la frontière 


système-milieu extérieur. 


Volume indéformable 


Cf 


—« Flux de 


Flux de o o | matière 
matière "= UEc, Ep mm M; 


mi (R -KU : Ci) 






(Pi-V;Ui Ci) LY Lo 





Figure 2.7 : Schéma d’un système ouvert à plusieurs entrées (i) et sorties (j) 


Le bilan énergétique d’un système ouvert s’écrit donc : 


AE, = AU + AE + AEp — 


C2 C2 
entrée sortie 


Ou, W; est le travail technique (utilisable techniquement) 


Si on introduit l’enthalpie H=U+PV: 


C2 C2 
entrée sortie 


S1 de plus on note par : 


2 
hr = h+ — +gz:  l’enthalpie totale de l’unité de masse de fluide, [kJ/kg] 
La relation (2.3) devient : 


4 i ; ; ; ae ag ; r 
Le travail de transvasement est le travail requis pour faire passer la matière à travers les sections d’ entrées et 
de sorties du systeme ouvert (voir démonstration dans annexe) 
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AE, = Q + W, + > T ee > m, (hr), (2.3a) 


entrée sortie 


ou encore sous forme différentielle : 


dE, = 60 + 6W, + > dm; hr, — > dmj hr (2.3b) 


entrée sortie 


Ces équations sont valables quel que soit le nombre de flux de matière (c'est-à-dire le nombre 
des entrées et sorties), que le régime soit permanent ou transitoire. 
En divisant (2.3b) par le temps dt, on obtient une écriture en puissance du bilan d’énergie: 


Èr = Q+W, + > m hr, — 2 mj hr; 


entrée sortie 
avec 
Èr la variation de l’ énergie totale du système pendant l’ unité de temps, [kW] 
Q la puissance thermique échangée entre le système et le milieu extérieur, [kW] 
W, la puissance technique échangée entre le système et le milieu extérieur, [kW] 


m;,m, le débit massique du fluide, [kg/s] 


Dans de nombreuses applications, le système ouvert possède une seule entrée (1) et sortie (2) 
(Fig. 3.8). On a alors, en considérant le principe de conservation de la masse m, = m, =m : 


Ep = Q + W; + m(hr, — hr,) 


Volume indéformable 








Z mı 
Système ouvert | (2, P, U5, C3 
U, Eç, Ep ds 
m 
(Pi. Vi Un Cù 4 


as 


Figure 2.8 : Schéma d’un systeme ouvert a une seule entrée et sortie 


Si le système ouvert est en régime permanent (pas de variation des propriétés au cours du 


temps), son énergie totale est constante : 
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Èr = dEr = 0 
Le premier principe de la thermodynamique (en terme de puissance [kW]) se réduit donc a: 


Q + W, = (hr, SE hr.) = mAhr 


2 2 


, Cz — Ci 
=m (ha — hy) + — > — t g(Z2 — 21) 


en divisant cette relation par le débit massique m, on aura : 


q +w, = Ahr; , [kJ/kg] 


X\ 


w+ et q sont, respectivement, les échanges d’énergie mécanique et thermique rapportés à 
l’unité de masse qui entre et sort du système. 


En multipliant par la masse m, on obtient : 


Q +W, = AH; , [kJ] c 


ou encore pour une transformation élémentaire : 


ôQ + OW, = dHy 


2.4.1 Application du premier principe à quelque machines en régime permanent 

> Compresseur et Pompe 

Un compresseur (Fig. 2.9) est un organe qui transfère de l’énergie (travail technique W,) à un 
gaz, ce qui augmente la pression de ce dernier (P2 > Pı). Une pompe fonctionne aussi de la 
même façon mais en utilisant un liquide. Si on suppose négligeable (et c’est pratiquement le 
cas) l’échange de chaleur avec le milieu extérieur (Q = 0), la variation de l’énergie cinétique 


Cs" 


Be 
(AE =m ana = 0) et la variation de l’énergie potentielle (AE, = mg(z> — Z1) = 0), le 


bilan d’énergie s’écrit des façons suivantes : 


W, = AH = H, — H, > 0 [kJ] 


W, = mAh = m (h, -h )>0 [kW] 
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Figure 2.9 : Schéma d’un compresseur 


Le travail consommé par le compresseur (ou la pompe) se retrouve dans l'augmentation de 


l’enthalpie du fluide (Hz > H;) 


> Turbine à gaz (ou à vapeur) 
Dans cet organe (Fig. 2.10) le fluide (gaz ou vapeur d’eau à haute température) fournit de 
l’énergie (W,) au milieu extérieur. Il en résulte une chute de pression du fluide entre l’entrée et 


la sortie de la turbine (détente). En négligeant Q, AE-et AE, on aura: 


W, = mAh = m (h, —h,) < 0 [kW] 





Figure 2.10 : Schéma d’une turbine 


Le travail fournit au milieu extérieur provient de la chute d’enthalpie du fluide (Hz < H;) 


> Echangeur de chaleur 


Un échangeur de chaleur (Fig. 2.11) transfère de l’énergie thermique (Q) d’un fluide chaud 
(A) a un fluide froid (B). Le système composé des deux fluides (A et B) est supposé 
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adiabatique et la variation de l’énergie cinétique et potentielle est nulle. De plus, comme il n’y 
a pas de mise en jeu de travail technique dans ce type d’organe, le bilan d’énergie (pour tout 


l’échangeur) s’écrit simplement: 
0 = m,Ah, + mgAhz 


[kW] 
= ma (haz — Nar) + Mp (hg — hga) 





ts Ls 


Figure 2.11 : Schéma d’un échangeur de chaleur 
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Chapitre 3 : Les gaz parfaits 


3.1 Définition d’une équation d’état 

Une équation d’état est une relation de dépendance entre les différents paramètres d’un 
système en état d'équilibre. En thermodynamique, elle s’écrit sous la forme générale 
suivante : 

f(P,V,T,n) =0 

L’équation d’état peut être appliquée : 

- le long de la transformation si celle-ci est réversible ; 

- seulement pour l’état initial et final (a condition qu'ils soient des états d’équilibres) si 


la transformation est irréversible. 


3.2 Gaz parfaits 
Gaz hypothétiques dont les particules sont considérées comme des points matériels! sans 
interaction (attraction ou répulsion) entre eux. Les conditions thermodynamiques pour 
lesquelles un gaz quelconque se comporte comme un gaz parfait sont les suivantes : 
- très faible pression (P > 0 bar) 
- faible densité (gaz dilué) 
- faible pression (P << Per) et haute température (T >> Ta 
Par exemple à T = 20°C et P = 1bar : 
o l’Hydrogène (Te = -240°C, Per = 13bar), Oxygène (Te = -119.2°C, 
Pa = 50.5bar) et l’ Azote (Ter = -147.2°C, Per = 33.9 bar) se comportent comme 
des gaz parfaits 
o Par contre, les gaz lourds ou denses comme la vapeur d’eau (Ter = 374.2°C, 
Pr = 220.9bar) ou les réfrigérants vapeurs (exp: R134a; Te = 100.9°C, 


Pa = 40.7bar) ne sont pas considérés comme des gaz parfaits. 


' Cela veut dire aussi que le volume propre des molécules d’un gaz parfait est infiniment petit par rapport a 
l’espace dans lequel elles se déplacent (c'est-à-dire par rapport au volume occupé par le gaz). 
* Ces deux conditions doivent être vérifiées à la fois 
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3.3 Lois essentielles des gaz parfaits 
3.3.1 Loi de BOYLE-MARIOTTE : (loi de compressibilité des gaz) 

Boyle et Mariotte établissent qu’à température constante, le produit de la pression (P) par 
le volume (V) d’une masse de gaz invariable est une constant. : 

P xV = Cst (3.1) 

Cette relation constitue l’ équation caractéristique d’une transformation isotherme d’un gaz 
parfait. 
Expérience : considérons un mécanisme cylindre-piston renfermant un gaz parfait maintenu à 
une température constante. En faisant varier la pression, le volume varie dans le sens inverse 
(P est inversement proportionnelle à V). Sur le diagramme de Clapeyron (P-V), la loi de 


variation de P en fonction de V constitue une branche hyperbolique (y = a/x)(Fig.3.1) 


Pa 
_ Cste 
y 





PCV) 


— Figure 3.1 : Loi de Boyle-Mariotte 
— ; (P x 1/V) 


Vl V2 V3 y 











En tout point de la transformation réversible représentée par la Figure 3.1 la relation (3.1) est 
vérifiće. Nous pouvons donc écrire : 


Pi XV, =P, xV, =- =P, XV, = Cst 


3.3.2 Loide GUY LUSSAC : (dilatation des gaz à pression constante) 
A pression constante, le volume d’une masse invariable de gaz est proportionnel à sa 
température absolue : 
V = Cst XT ou V/T = Cst (3.2) 
Cette relation est l’ équation caractéristique d’une transformation isobare d’un gaz parfait. 
Expérience : considérons le même mécanisme cylindre-piston. La pression est maintenue 
constante sur le piston. En faisant varier la température du gaz, le volume varie 


proportionnellement. La loi de variation de V en fonction de T est linéaire (Fig.3.2). 
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Figure 3.2 : Loi de Guy Lussac (V « T) 





Nous pouvons donc écrire : 

Vu Fy V, 

Lea a a 

L h Ty 
3.3.3 Loi de CHARLES AMONTONS (variation de la pression à volume constant) 
A volume constant, la pression d’une masse invariable de gaz est proportionnelle a sa 
température absolue. 

P=Cst XT ou P/T =Cst (3.3) 

Cette relation est l’équation caractéristique d’une transformation isochore d’un gaz parfait. 
Expérience : Considérons un réservoir fermé, à parois fixes, rempli d’un gaz parfait. En 


faisant varier la température, la pression varie proportionnellement. La loi de variation de P 


en fonction de T est linéaire (Fig.3.3). 





isochore ¢ 
































Tl T T3 + Figure 3.3 : Loi de Charles (P « T) 


Là aussi, nous pouvons écrire : 


P P P 
res eee ee Cst 
ae T, 
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3.3.4 Loi d’'AVOGADRO 


Avogadro établit que dans les mêmes conditions de température et de pression, le volume 
occupé par une mole d’un gaz parfait (volume molaire v = V /n)) est le même quelque soit le 
gaz considéré. Autrement dit à 7 et P constants, le volume (V) est directement proportionnel 


au nombre de moles (n) : 


V=Cstxn (3.4) 
Ou encore 
V. V. V 
duo 
M MM Nn 


v AT=0°C (273,15 K) et P= 101,3 kPa (condition normale) le volume molaire de tous 


les gaz est V = 22,4 m°/kmol 


3.4 Equation d’état des gaz parfaits 
La combinaison des lois précédentes donne lieu a une relation plus générale dite équation 
d’état des gaz parfaits et qui s’écrit : 
f(P,V,T,n) =PV-nRT=0 
ou tout simplement : 
PV=nRT (3.5) 
avec 
P : la pression absolue, [kPa] ; 
V : le volume occupé par le gaz, [m°] ; 
n : le nombre de moles, [kmol] ; 
R : la constante universelle des gaz parfaits, =8.314 kJ/kmol.K ; 


T : la température absolue, [K]. 
L’équation d’état (3.5) peut s’écrire de différentes façons : 


v divisons la relation (3.5) par n, pour avoir : 
Pov =RT (3.6) 
v : le volume molaire = V/n, [m°/kmol] 
v en écrivant dans (3.5) que n = m/M, on obtient : 
PV=mrT (3.7) 
M : la masse moléculaire du gaz, [kg/kmol] 


m : la masse du gaz, [kg] 
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r : la constante du gaz parfait = R/M , [kJ/kg.K] 
v divisons (3.7) par la masse V : 
P=prT (3.8) 
p : la masse volumique = m/V, [kg/m°] 
v divisons (3.7) par la masse m : 
Prv=rT (3.9) 


v : le volume massique = 1/p, [m°/kg] 


3.5 Lois de Joule 
Expérience de Joule : Joule dans son expérience (Fig. 3.4) considère deux réservoirs 
(plongés dans un bain d’eau a To) reliés entre eux par un robinet fermé. Dans le 
réservoir (A) il place une quantité d’air comprimé a Po (gaz parfait) et dans l’autre (B) 
il crée le vide (P = Obar). Lorsqu'il ouvre le robinet, Il’ air passe dans (B) et se détend : la 


pression de I’ air à diminuée (P;< Po) et le volume qu’il occupe a augmenté (V;> Vo). 


|| Thermomètre || Thermomètre 





AIR reservoir 








Figure 3.4 : Expérience de Joule 


Joule constate que la température To du bain reste constante le long de l’expérience : il n’y 
a donc pas de chaleur échangée entre lair et l’eau (Q = 0). Il n’y a pas non plus de travail 
mis en jeu puisque rien n’est mis en mouvement, ni soulevé, ni entrainé..... et donc dW = 0 
aussi. Le premier principe nous permet d’écrire : 
dU = óW + ôQ = 0 
Joule conclut que l’énergie interne (d’un gaz parfait) est constante à T constante même si la 


pression P et le volume V varient, d’où la première loi de Joule 
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e Première loi de Joule: L'énergie interne d’un gaz parfait ne dépend que de la 


- 3 
température . 


ur) > du = (9%) ar + (224 av = (2%) ar 
OT /y V/ OT/y 

0 
u est l’énergie interne massique [kJ/kg] 


(du/OT)y , on aura: 


comme la chaleur spécifique massique a volume constant cy 
(2) 
du=cydT > Au= | Cy AT = cy AT 
(1) 
v Variation de l’énergie interne massique u [kJ/kg] 
Au = Uz — Uz = Cy AT 
v Variation de l’énergie interne molaire u [kJ/kmol] 
Au = Uz — U = Cy AT 
Cy est la chaleur spécifique molaire a volume constant [kJ/kmol.K] 
v Variation de l’énergie interne U [kJ] 
AU = U, —U, =m cy AT = nc, AT 


m et n désignent respectivement la masse et le nombre de moles du gaz 


e Deuxiéme loi de Joule : L’enthalpie d’un gaz parfait ne dépend que de la température. 
Démonstration : considérons l’enthalpie massique A [kJ/kg] 


Puisque h=u+Pv etque Pv=rT => h=u(T)+rT = A(T) 


hr) = an = (Z) ar + (Ê ) dP = (=) aT 
~ \OT/ > P/r  \OT/p 


0 
comme la chaleur spécifique massique à pression constante cp = (0h/0T) , on aura: 
(2) 
dh=cpdT > Ah= | Cp dT = cp AT 
(1) 
v Variation de l’enthalpie massique h [kJ/kg] 
Ah=h;-h = cp AT 


v Variation de l’enthalpie molaire h [kJ/kmol] 


? Ce qui veut dire que l’énergie interne d’un gaz parfait est uniquement composée de l’énergie cinétique 
microscopique des molécules, liée à l’agitation thermique (l’énergie potentielle microscopique d’interaction est 
donc exclue). 
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Ah = h, — h; _ Cp AT 
Cp est la chaleur spécifique molaire a pression constante [kJ/kmol.K] 


v Variation de l’enthalpie H [kJ] 
AH = H, — Hı = m Cp AT = Nn Cp AT 


3.6 Transformations réversibles d’un gaz parfait 

On s’intéressera dans ce paragraphe à la représentation, sur le diagramme de Clapeyron, 
des transformations thermodynamiques des gaz parfaits, dans le cas d’un système fermé et 
ouvert. On établira également les expressions du travail des forces de pression, du travail 


technique et de la chaleur mises en jeu par les systèmes en question. 


3.6.1 Transformation isotherme (7=cste) 

Une transformation isotherme est difficile à réaliser sauf si on place notre système dans un 
bain thermostaté (à 7=Cste) tel que la glace fondante ou l’eau bouillante... 
L’équation caractéristique d’une telle transformation, lorsque le gaz est parfait, est donnée par 
la loi de Boyle-Mariotte vue précédemment : 


PV =cste = nRT 
a/ Systeme fermé 


Travail des forces de pression (W): 
(2) (2) 
nRT V> 
W = — f PdV = — f (F) av = -nrTin (2) 
V Vi 
(1) (1) 
pour le travail spécifique massique (w = W/m): 


= rtn (7) 
W = -r AVA 


Figure 3.5 : Travail de forces de pression 
d’une compression isotherme 
dans un système fermé 





Quantité de chaleur (Q) : 
comme le gaz est parfait AU (T) et isotherme (AT = 0), alors: AU = 0 = 
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V. 
Q+W=AU=0 > Q = -W = nRTIn (32) 
1 


ou encore pour une quantité de chaleur rapportée à l’ unité de masse 


= rrin (7) 
q=1Tin\7 


1 


b/ Système ouvert 


Travail technique (W): 
(2) (2) 
nRT P; 
W, = f VdP = f (=) dP = nRTin (=) 
P P, 
(1) (1) 


rapporté à l’unité de masse 


P, 
w+ = rTin (=) 
comme la transformation est isotherme 
PV, = PV => ae 
P Vz 


en remplaçant cette relation dans l’expression de Wt et w;, on obtient : 


V2 
W; = —nRTIn (2) = W 
vi 


et Wt = W 






T= este 


Figure 3.6 : Travail technique d’une compression 
isotherme dans un système ouvert 


Quantité de chaleur (Q) : 
à partir du 1°” principe de la thermodynamique pour un système fermé (voir CH 2): 
Q = —W, + AH, 


ou encore q = —w, + Ahr 


3.6.2 Transformation isobare (P=cste) 


L’équation caractéristique d’une transformation isobare d’un gaz parfait, est donnée par la 


loi de Gay Lussac : 
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a/ Système fermé 


Travail des forces de pression (W): 
(2) 
W = — f PdV = —PAV = —P(V, — V1) 
(1) 
ou w = —P (v, — v) 


v est le volume massique [m°/kg] 


l P=cste 2 


Figure 3.7 : Travail de forces de pression 
d’un échauffement isobare 
dans un système fermé 





Quantité de chaleur (Q) : 
— (U> + PV,) — (U, + PV;) 
= H, — Hı = AH = mcpAT 
ou encore pour une quantité de chaleur rapportée à l’unité de masse 
q — Ah = CpAT 
b/ Système ouvert 


Travail technique (W): 


Figure 3.8 : Travail technique d’un échauffement 
isobare dans un système ouvert 
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Quantité de chaleur (Q) : 


Q=AH, et q= Ahy 


3.6.3 Transformation isochore (V=cste) 
L’équation caractéristique d’une transformation isochore d’un gaz parfait est donnée par la 


loi de Charles Amontons : 


a/ Systeme fermé 


Travail des forces de pression (W): 
(2) 


w=- | Pav =0 d wai 
(1) 


Figure 3.9 : Travail de forces de pression 
d’un échauffement isochore 
dans un systeme fermé 





Quantité de chaleur (Q) 
Q=AU=mcyAT et q = Au = cyAT 


b/ Système ouvert 


Travail technique (W): 


(2) 
W; = | VdP = VAP et Wt = vAP 
(1) 


Figure 3.10 : Travail technique d’un échauffement 
isochore dans un système ouvert 
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Quantité de chaleur (Q) : 


Q = —-W, + AH; et q = —w, + Ahr 


3.6.4 Transformation adiabatique’ (O=0) 
Cherchons |’ équation caractéristique d’une transformation réversible d’un gaz parfait : 
le premier principe —> dU = ôW (a) 
gaz parfait > dU = ncydT (b) 
trans.réversible > W = —PdV = —(nRT/V)dV (c) 


en remplaçant (b) et (c) dans (a), on aura : 


dT RAV _ DY dr. po 
r> gy YOY 7t+aY-)7= 


y est la constante adiabatique (= cp/Cy = Cp/Cy) 
en intégrant cette relation, on obtient : 
nT +InVY-Y) =cste = 

TVY-VD = cste > T0 D = TV, — 1° relation de Laplace 
En utilisant l’équation des gaz parfaits pour éliminer 7 de cette relation, on obtient la 
deuxième forme de l’équation caractéristique d’une transformation adiabatique réversible : 

PVY =cste = P,V," = PV” — Loi de Poisson 

Faisons la même chose avec V, on aura: 


1-y 


rpl y Ve cste => rp Cr) = T p, 7) — 2®™e relation de Laplace 
a/ Système fermé 


Travail des forces de pression (W): 
(2) (2) 


cste cste 
-e 2 _ a-y) _ ya”) 
W = | Pav = [jw T lh yo | 
(1) (1) 


comme la constante cste = P,V,’ = P>V," , alors 





W = ELUS pb pr) PV aw = 1 
(y—1)l* : œ —-1[\" 


4 


Une transformation adiabatique et réversible est isentropique (à S=cste) 
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| PV, JT. 
1°" relation de Laplace — — a | 
v- DIT 
1-y 
2° relation de Laplace — W = (=) —1 
 —1)|\P 


On peut également retrouver l’expression du travail a partir du premier princioe de la 


thermodynamique, tel que : 


W = AU = mc,AT = ncy AT 





Su = Fi tate gl 

sachant que : Cy = = et Cy= PE 

W å AT AT 
= M — = 7) 
Yt) (y — 1) 
b/ Système ouvert 
Travail technique (W): 
dP dV 
PVY =cste => InP +ylnV = cste => ty 0 => 


VdP = —-yPdV => ÔW = yoW = W, = yW 


Donc 
1-y T. 
= yY V2 W, = —— pv 2-1 
= pen rv. (E) 1 P-L, 
p 1-y 
Y 1\ Y 
jeter (=) si 
í (y — 1) a P, 


Le travail technique peut être calculé à partir du 1° principe de la thermodynamique, tel 
que p p p p y q 


que : 
W; = AHr = AH + AEp + AEc 
si AEp = 0 et Eç = 0, alors : 


W, = AH 


PO 7 


Ce résultat a déja été retrouvé au chapitre précédent 


5 ai à : C Cp Multi ty. ee : 
Si on combine ces relations : F = = = y (définition) etcp —cy =r , Cp — Cy =R, on aura facilement: 
V 
be oP. oe 
Vi y-1 V y-1 
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La pente de l’adiabatique et de l’isotherme dans le digramme (P-V) ? 


La pente d’une quelconque courbe dans le diagramme (P-V) est donnée par dP/dV 


L’adiabatique L’isotherme 
PVY = cste => InP + yinV = cste = PV =cste = lnP + InV = cste > 
salen a _, oP . 2s dP dV dP P 
pity” dv. vy pty" y 
hs P P 
pente de l'adiabatique = -Y7 pente de l'isotherme = 7, 


pente adiabatique 
—— =y>l 
pente isotherme 


On conclut que dans le diagramme (P-V), l’adiabatique est plus raide que l’isotherme 


(Fig. 3.11). 






__ isentropique k = Y 
\— isotherme k = | 


Figure 3.11 : La pente de l’adiabatique 
D — plus importante que celle 
de l’isotherme 


3.6.5 Transformation polytropique 
La transformation polytropique s'approche davantage d'une transformation réelle, son 


équation caractéristique est : 


1-k 
PVk =cste , TV“’-Y) =cste , rp) = cste 


k est le coefficient polytropique. En faisant varier ce coefficient de 0 — œ on décrit toutes 


les transformations possibles (Fig.3.12) : 


nn 
op me ue 


DCS PV® = cste = P@/V = cste isochore 
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Figure 3.12 : Transformations polytropiques 
k= dans le diagramme de Clapeyron 





a/ Systeme fermé 


Travail des forces de pression (W): 


Reprendre les expressions du travail pour une transformation adiabatique et remplacer y 














par k : 
Ps [MY PV D 
-&-5[(w) j a 1 
P, V- P. = W = j AT = R A 
aea rairai 


Quantité de chaleur (Q) : 


R 
ôQ = dU — dW = ncydT —n——~dT 


(k — 1) 
comme Cy = =n alors : 
Yo (y-1) 
(k — y) _ 
6Q = nR——_—_dT = nc, dT 
(k—1)(7 - 1) i 
avec : 
: (K-y) _(k-y) We 
k RDO- D) K-D” 
d’où 
Q = nc, AT = mc, AT 
avec : 


— (K-y) _(k-y) Y) 
kl RDO- D) &-D” 
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expression de 59 
isochoe (k=) | ay | OQ = mgd = dU 


isobare (k = 0) = = ôQ =nCpdaT = dH 


adiabatique (k = y) 


b/ Systeme ouvert 





Travail technique (W): 
Reprendre les expressions du travail technique pour une transformation adiabatique et 


remplacer y par k : 











W; = kW 
1—k PV, [IT 
W. =k Pi Vi 5 24 W, = 1V1 2-1 
t "(k= DIV, (k — 1) IT; 
1—-k kr kR 
P,V. Py\ k = m—— AT = n———- AT 
We = kay (=) a LL E “= 
Quantité de chaleur (Q) : 
Q = —W, + AH; 
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Tableau récapitulatif des principales énergies liées aux transformations réversibles 
de gaz parfaits dans un système fermé et ouvert 


Transformations isotherme isochore isobare adiabatiqu polytropique 
réversibles (T=cste) (V=cste) (P=cste) (O=0) 


e 
PVY = cst PV* = cst 
Equation PV = cst P/T = cst V/T = cst 
d / / TVY-) = cst 
TPO) = cst mt al k)/k — cst 


nR 
—nRTIn (— —PAV AT = A 
eat ED 
fermé 





TV ED = cst 











mc, AT = N Ck AT 
t 


Af = fs — fı (1: état initial, 2 : état final) ; AH; =AH+AEc+AEp ; = De, ; a=% ; R= 8.314 kj/kmol K 
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Chapitre 4 : Mélange de gaz parfaits 


4.1 Mélange idéal de gaz parfaits 

Un mélange de plusieurs gaz parfaits est dit idéal s’il n’y a pas d’interaction physique ou 
chimique entre les différentes molécules de gaz mélangés. Le mélange possède alors les 
propriétés suivantes : 
- son énergie interne (U) est la somme des énergies internes (Ui) des gaz qui le composent. 


Il en est de même pour son enthalpie (H) et son entropie (S) : 


U= NU, H=) H; #35 
i=1 i=1 i=1 


- il obéit à la loi de joule : U=U(T) et H = H(T) 


- son équation d’état est : PV =nkRT 
4.2 Composition d’un mélange de gaz parfaits 


Soit un mélange, dans l’état (P,V T), de k gaz parfaits. Chacun de ces gaz, repéré par un 


indice i dans le mélange, est présent à raison de n; moles (Fig. 4.1). 


Amoles 
mélange 





| 
l 
| 
| 
| 
| 
l 
| 
| 
| 
l 
| 
_ 


Figure 4.1 : Mélange de plusieurs gaz parfaits 


Volume = V 


Le nombre total de moles dans le mélange est la somme des moles de chacun des gaz qui le 


compose : 
k 
n = » nN; 


i=1 
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De la même façon, la masse totale du mélange est la somme des masses de chacun des gaz qui 
le compose : 


i=1 
Pour connaitre la composition de ce mélange, deux grandeurs sont généralement utilisée : la 


fraction molaire ou la fraction massique, définie comme suit : 


= La fraction molaire (x;) du gaz i dans le mélange est égale au rapport du nombre de mole 


n; de ce composant sur le nombre total de mole du mélange : 


k 
x ey 
nr avec x; = 1 


i=1 


N; : le nombre de mole du gaz i, [kmol] 
n : le nombre total de moles, [kmol] 


= La fraction massique (y;) du gaz i est égale au rapport de la masse m; de ce composant 


sur la masse totale du mélange : 
k 
M; 
re <1 avec tHe 
i=1 


m, : la masse du gaz i, [kg] 


m : la masse totale du mélange, [kg] 


On peut établir une relation entre la fraction massique et la fraction molaire, telle que : 


Cm nM, nM; — Mi 
om nM nM UM 





avec : 
M; : la masse molaire du gaz i, [kg/kmol] 


M: la masse molaire moyenne (ou apparente) du mélange [kg/kmol], dont le 


calcul sera présenté ultérieurement. 


4.3 Pression partielle - volume partiel 
= La pression partielle P; du gaz i dans le mélange est la pression qu’exercerait ce gaz s’il 
était seul a occuper, a la température 7, le volume V offert au mélange (Fig. 6.2a). Le gaz i 


étant parfait, ona: 
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P; V = n;RT 
étant donné que le mélange est idéal : 
PV = nRT 
on en déduit : 
Ni 
P; = ?— X;P 
n 


on peut tout à fait vérifier que : 


> =P 


i 
C’est la loi de Dalton : La pression P d’un mélange idéal de gaz parfaits est égale à la 


somme des pressions partielles des différents constituants. 


Le volume partiel V; du gaz i est le volume qu’occuperait ce gaz à la température T et 


pression P du mélange (Fig. 6.2b) : 


PV, = n;RT (7) 
or le mélange est idéal : 
PV = nRT 
on en déduit que : 
V; = y = x;V 
n 


et 


Mav 
i 


C’est la loi d’Amagat: Le volume V d’un mélange idéal de gaz parfaits est égal à la 


somme des volumes partiels des différents constituants. 


. Température du a Température du 
IN mélange (7) Pression du 
_ mélange ( P) 






mélange (J) Pression partielle (FP; ) 


| 
CET << 2. 2. 2 il 








nn uen 
Volume du 


(a) mélange (F) 


Figure 6.2 : Schéma montrant (a) la pression partielle et (b) le volume partiel 
d’un composanti du mélange 
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44 Calcul des propriétés d’un mélange idéal de gaz parfaits 
On se propose de déterminer les propriétés d’un mélange de gaz parfaits, connaissant sa 


composition et les propriétés des gaz parfaits qui le constituent. 


= Masse molaire moyenne (M) 


m : M; n;M; 
re dim | LimiMi =\ xm, 


n n n 





i 
Exemple : supposons que la composition molaire de l’air est de 79% de N2 et 21% deO». 
calculons sa masse molaire a partir des propriétés thermodynamiques de ses composants : 


My2=28 kg/kmol Mo2= 32 kg/kmol. 


Solution : 
Mair = X02Mo2 + Xn2Mnz = 0,21 x 32 + 0,79 x 28 = 28,84 kg/kmol 


Remarque : Dans la littérature Mair = 28,97 kg/kmol. La différence entre cette valeur et 
celle qu’on vient de calculer est due au fait d’avoir négligé les autres composants de I’ air (tel 


que CO: ....) 


= Energie interne, enthalpie et entropie massiques (u,h,s) 


U= dus mu= > mu, su=) yu 
i i i 
H=) H> mh= > mh, sh=) yh 
i i i 
s=% S> ms =% ms >s= Y yisi 
i i i 


Ui, hi et s; désignent respectivement l'énergie interne massique, l’enthalpie massique et 


l’entropie massique du gaz parfait i. 


= Chaleur spécifique massique à pression ou volume constant (Cp, Cy) 


dH; 
l 
Cp LS, CPi 
i 


dU; 
l 
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(y= > Yi Cvi 
i 
= Chaleur spécifique molaire à pression ou volume constant (Cp, Cy) 
dH; _ Des 
H = DH =o - ar UPS QUE © 


l 
Cp = Xi Cp; 


i 


dU; E _ 
U = 25 =o = AT = Du, > 
L 
Cy = x, 


i 








= Constante adiabatique y 
La constante adiabatique du mélange est donnée par le rapport des chaleurs spécifiques à 


pression et volume constants, calculées ci-dessus : 
r- Cp _ Cp 
Cy Cy 


Po 7 


ou encore a partir de la relation suivante (déja établie) : 


Gy =) xt, > = x — 


L 


: +1 
= 7 
2 
Ll 
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Chapitre 5 : Les gaz réels 


5.1 Introduction 

L’équation d’état des gaz parfaits est un modèle mathématique très simplifié qui ne peut 
reproduire le comportement des gaz que dans un domaine de pression et de température très 
limités. Depuis, plus de 150 équations d’état (de gaz réels) empiriques ont été proposées, sans 
qu'aucune d’entre elles ne soit à la fois suffisamment précise et suffisamment maniable. La 
plus connue et la plus simple des équations d’état des gaz réels est celle de Van Der Waals 
(1873). Elle ne concorde pas de façon parfaite avec les résultats expérimentaux, mais toutefois 
elle permet d’analyser les phénomènes de changement de phase des gaz réels. Nous pouvons 
citer d’autres équations plus complexes telles que : l’équation de Readlich-Kwang, de Beattie 


Brigeman, de Benedict Webb Rubin..... 


5.2 Equation d’état de Van Der Waals 

Lorsque nous avons considéré les gaz parfaits, le volume des molécules ainsi que les 
phénomènes d’interactions existants entre elles ont été négligés. Van der Waals dans son 
modele (plus performant que celui des gaz parfaits, mais pas le meilleur) a tenu compte du 
volume des molécules et de leurs interactions mutuelles. Les conséquences directes de ces 


hypotheses sur la relation PV = nRT sont les suivantes : 


e 1” conséquence 
Plaçons n molécules d’un gaz dans un récipient de volume V. Chaque molécule 
possède un volume propre non nul b, appelé covolume. La quantité nb représente alors 
le volume inaccessible aux molécules. Autrement dit, le volume libre (disponible) au 
déplacement des n molécules n’est plus V mais (V-nb). L’équation d’état devient 
(provisoirement) : 
P(V — nb) = nRT (5.1) 
2°" conséquence 
Comme chaque particule de gaz est attirée par les autres particules qui l’entourent, 
nous distinguons deux situations (Fig.5.1) : 
- toute particule loin des parois est attirée dans toutes les directions de la même 
façon (hypothèse d’homogénéité) et la résultante de ces forces d’attraction est 


nulle 
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- par contre, toute particule à proximité des parois subit une attraction 
dissymétrique (dirigée vers l’intérieur du gaz) qui réduit (par rapport au cas du 
gaz parfait) la pression contre les parois d’une valeur de a(n/V)? (appelée 
aussi pression interne) : Pep = Per + a(n/V)* 


Ainsi nous obtenons 1’ équation dite de Van Der Waals qui s’écrit! : 


r 
[> + az) (V — nb) = nRT (5.2) 
Ou aussi sous la forme : 
pe nRT n? 4 
~(V—nb) “V? ped) 


avec 


P : la pression absolue du fluide (gaz réel) ; 

T : la température absolue ; 

a et b: des constantes positives de l’équation de Van Der Waals qui 
dépendent de la nature du fluide (et non pas des variables d’état P,V....). Le 


tableau 1 donne quelques valeurs de ces paramètres. 


Tableau 1 : Constantes a et b de l’équation de Van Der Waals 


Fluide a x 10 b x 10° 
[Pa.m°/mol”] [m°/mol] 
H20 5,536 3,049 





Figure 5.1 : Illustration de l’effet de l’interaction moléculaire 
sur la pression du gaz réel 


* Lorsqu’on spécifie P et T, l'équation de Vander waals est cubique en volume 
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Voici d’autres formes de l’équation de Van Der Waals : 


"4 (P+a 3) v - nb) =nRT => P+ (—-b)n=ner = 
n 
n2 
(P+= —) (6 —b) = RT (5.3) 


RT a 
ou encore PE Gp) = 72 (5.3a) 


v : le volume molaire (=V/n) 
m\2 

P _ (ii) 
P+ar (V—nb) =nRT =|P+a y2 





(p+ CEE (5.4) 
TT a’ 
— (5.4a) 


Sencar. = (v—b') v? 


v : le volume massique (=V/m) , 
les coefficients de l’équation de Van Der Waals, tels que a’ = a/M* et b' = b/M 


a'et b': 
Dans le paragraphe 5.4 nous verrons comment évaluer les paramètres a’ et b’(ou a et b) a 


partir des variables d’état (pression, température et volume) du point critique 


5.3 D’autres équations d’état de gaz réels 
Equation de Redlich-Kwong (à 2 coefficients) 








= RT a 
(5-b) 5(5+b)T/2 
2m5/2 
= 0,42748 —— t b= 0,08664 
a po P. 


C 


Te et Pe sont la température et la pression au point critique 


Equation de Beattie-Bridgeman (à 5 coefficients) 
RT E y ô 
P=—+—+—+— 
DUO T VV. 


avec 
B = BRT — A — cR/T? 
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y = —BbRT + Aa — BcR/T? 
ô = BbcR/T* 


A, B ‚a, b etc sont des coefficients (tabulés) de l’équation 


e Equation de Benedict Webb Rubin (a 8 coefficients) 


RT C\ 1 (bRT-a) aa C y 5 
— E ae E ee 8 ee Pe |! \o-y/v 
P A (BRT A =) = 7 oes (1 4 zz) 5 


A, B, C, a, b, c, a et y : les coefficients (tabulés) de l’équation 


5.4 Isothermes de Van Der waals 
La représentation des isothermes de Van Der Waals sur le diagramme de Clapeyron (P,v), 
nous permet de constater les faits suivants (Fig. 5.2): 

v Pour des températures constantes élevées, P et v sont grands par rapport aux termes 
correctifs a/v* et b = nous retrouvons sensiblement Pv = rT (comportement d’un 
gaz parfait) 

vA mesure que les températures baissent les isothermes se déforment 

vV Lorsqu’on atteint une température dite critique (Tc), l’isotherme présente un point 
d’inflexion (et de plus à tangente horizontale) 

v En dessous de cette température (T< Tc), les isothermes présentent une zone 
d'instabilité mécanique qui traduit l’existence de deux phases fluide (liquide et 
vapeur). Mais expérimentalement, on observe des paliers à pression constantes 
(segment [ED] sur Fig.5.1) dits paliers de liquéfaction (ces cas seront traités 
ultérieurement dans le chapitre sur les changements de phases) 


Comme au point critique la tangente est horizontale et qu’il y a un point d’inflexion, nous 


avons la dérivée première et deuxième de P par rapport à v qui s’annulent: 


OP 
(—) = 0 (tangente horizontale) (5.5) 
OV/r=Te 

0°P ; jà 

TE =) (point d'inflexion) (5.6) 

T=T¢ 
Soit : 
rT, 2a’ rT 2a 


= 0 (5.7) 


bte Ub v 
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2rT. 6a’ n rT 3a’ (5.8) 
$$ c — —= => — = — . 
(we = b')? Ve“ (we = b')? v,4 
P., Vc, Te désignent respectivement la pression absolue, le volume massique et la température 


absolue au point critique. 


En divisant membre a membre les égalités (5.7) et (5.8) nous obtenons : 


Ve = 3b' (5.9) 
La substitution de (5.9) dans (5.7) puis dans (5.4a), donne respectivement: 
T= LS (5.10) 
27b'r 
a! 
e— 27p2 (5.11) 


Pression (atm) 





300 
Ving fom? . mot") Figure 5.1 : Isothermes de Van Der 


— otherness expeérimentaies 


— — isotherme calculée d'après la loi des gaz partaits waals dans le diagramme de 
-= sotherme calculée d'après l'équation de Van cher Waals 


Clapeyron 


Donc réciproquement, à partir de la mesure de (P., v., Te) nous pouvons calculer les 


coefficients a’ et b’ (ou a et b) de l’équation de Van Der Waals, tels que : 


v 
a' =3P.v2 et b'= = (5.12) 


même la constante r du gaz parfait peut être exprimée en fonction de (P, Ve, Te) (rapport de 


(5.10) sur (5.11)): 
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_ BU 
=T, 





(5.13) 


Remarque : 

Le comportement des gaz réels dévie de celui des gaz parfaits au voisinage de la région de 
saturation et au point critique. Cette déviation peut être quantifiée par un facteur dit facteur de 
compressibilité noté Z (c’est aussi un facteur de correction): 

z-o PV PVP? 
rT nRT mrt RT 
e siZ=1 — gaz parfait (Fig.5.2) 
e siZ+ 1 — gaz réel (Fig.5.2) 


Figure 5.2 : Le facteur de compressibilité Z 


de l’hydrogène en fonction de la pression 





pour différentes températures constantes 
0 100 200 


p (atm) 


5.5 Loi des états correspondants 


S1 nous substituons a’, b’ et r dans la relation (5.4a) : 
a’ 
(p+ S)w-09 + à à 
V 


par leurs expressions (5.12) et (5.13), il apparait une forme particulière de l’équation de Van 


Der Waals qui ne dépend d’aucun caractère physique de la matière (équation généralisée): 


(P, + 5) (3v, — 1) = 8T, (5.14) 


P. = P/P. : la pression réduite 
T, = T/T, : la température réduite 


v, = v/v, : le volume réduit 
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Cette relation nous permet d’énoncer le principe des états correspondants : Si deux gaz de 
Van Der Waals ont deux variables d'états réduites égales, alors la troisième variable d’ état 


réduite est aussi égale. On dit que les deux gaz sont dans des états correspondants. 


Exemple : considérons de l’hélium” à la température de 10,4 K (soit 7r =2) et pression de 
2,3bar (soit Pr=1), et de l'hydrogène” à la température de 66,4 K (soit Tr =2) et pression de 


13bar (soit Pr=1) ces deux gaz occupent le même volume réduit 


5.6 Détente de Joule-Thomson 

La détente de Joule-Thomson est une détente isenthalpique (à H =cste). En pratique, elle 
est réalisée à l’aide d’une conduite calorifugée où le gaz, sous certaines conditions 
thermodynamiques, voit sa pression et sa température diminuées” à la traversé d’un obstacle 


(paroi poreuse, vanne, filtre, orifice...) (Fig. 5.3). 


(a) vanne 


no 


a aca Figure 5.3 : Schéma de quelques organes de 


détente de Joule-Thomson 


(c) tube capillaire 


Démonstration : Ecrivons le 1™ principe de la thermodynamique pour le fluide qui traverse 
la paroi poreuse (Fig.5.3) 
Q + W; = AH + AEc + AEp 


Comme Q, W,, AEc et AEp sont nuls alors : 


AH = 0 = H = cste 


paroi calorfugee 


P2172 





parol poreuse Figure 5.3 Détente de Joule Thomson 


* Point critique de He : T, =5,2K , P.=2,3bar, v, = 0,233 m3/kg 
? Point critique de H2: T,=33,2K ,P.=13bar , v, =0,52m°/kg 
* Lorsque le gaz est parfait la température reste constante lors d’une détente isenthalpique car H(T) 
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La détente de Joule Thomson est en effet une transformation isenthalpique. Elle est très 
utilisée pour la liquéfaction des gaz et dans plusieurs dispositifs de réfrigération. Toutefois, 
dans certains cas (conditions de 7 et P initiales), la température du fluide à la fin de la détente 
peut rester constante ou encore augmenter (situations indésirables). Le comportement de la 


température, lors de cette détente est décrit par le coefficient de joule-Thomson définit par : 


(=) 
JT = aP), 
Ce coefficient est une mesure de la variation de la température avec la pression durant un 


processus isenthalpique (c’est aussi la pente de la courbe h= cste dans le diagramme T-P). 
A partir de la figure 5.4, nous constatons que : 


<0 détente avec augmentation de la température (T, > T;) 
Ur À = 0 détente sans variation de la température (T, = T;) 
> 0 détente avec diminution de la température (T, <T:) 
La courbe d’inversion de l’effet Joule Thomson regroupe l’ensemble des points maximums 


des isenthalpiques, elle est donc donnée par la relation ur = 0. Il est tout a fait claire 


(Fig.5.4) qu’un effet de refroidissement par détente de Joule-Thomson ne peut avoir lieu que 


si le fluide est en dessous de la température d’inversion maximale notée Tj, : 


e S111>1;, alors Tz >T; (échauffement du fluide) 


e S17,< Ti alors Tz < T, (refroidissement du fluide) 


T temperature d'inversion 


Pa maximale 
(ae ral 
M 
ks 
\ 
\ 
A | 
Hrr > : Lp, < 0 
| h = const. 
| 
I 
! 
l Courbe d' inversion 
# 
Figure 5.4: Les lignes isenthalpiques dans le 


diagramme T-P 
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Chapitre 6 : 
Coefficients thermoélastiques et calorimétriques 


6.1 Les coefficients thermoélastiques (ou thermomécaniques) 


Ces coefficients permettent de prévoir le comportement thermoélastique’ d’un système 
(fluide dans notre cas). On les définis à partir de l’équation d’état du système f(P,V,T) = 0. 


Ces coefficients sont positifs et sont au nombre de trois : 


> Le coefficient de dilatation isobare œ qui mesure la variation du volume (V) liée à la 


variation de la température (T) à pression (P) constante (Fig.6.1): 


a= Gr), (6.1) 


OT 


Figure 6.1: Schéma montrant le 
sens physique du coefficient de 
dilatation isobare d’une substance 





100 kPa 100 kPa 100 kPa 
| kg | kg | kg 
Sustance a faible Sustance a forte 
valeur de Q valeur de Q 


> Le coefficient d’augmentation de pression isochore 64% qui mesure la variation de P en 


fonction T a V constant : 


_ 1 (2) a5 
> Le coefficient de compressibilité isotherme Xr (ou de dilatation isotherme) qui mesure la 


variation de V en fonction de P a T constante : 


=: sa 


Remarques : Le signe (-) dans cette relation permet de définir un coefficient positif car V 
diminue lorsque P augmente. Pour un système incompressible (pas de variation de V avec P), 


nous aurons : x7 = 0 


En tenant compte de l’identité de Reich (voir annexe) : 


1 w s ° id IN . . . 

C'est-à-dire la variation du volume du système avec T et P, ou encore la variation de sa pression P avec T. 
2 2 2 e o . 41, 2 œ 

Appelé également coefficient de compressibilité isochore 
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(Z) : (—) . (5) 4 
OV); OT) p OP/y 


Gi). ay -g 


et du fait que 


Nous pouvons établir une relation entre les différents coefficients thermoélastiques, telle que : 


æ = PP Xr (6.4) 
Dans le cas d’un gaz parfait pour lequel l’équation d’état f(P,V,T) = PV —nRT = 0: 


1 1 
a=p=z et AT D 


6.2 Les coefficients calorimétriques” (ou thermométriques) 

Considérons un fluide homogène en quantité bien définie (n = cst). L’échange de chaleur, 
de manière réversible, entre ce fluide et le milieu extérieur peut être exprimé sous différentes 
formes (selon le couple de variables indépendantes choisi : T,V ; T,P ou P,V) faisant appel a 


six coefficients calorimétriques généralement variables: Cy, Ch,l,h,À2,u 


> Variables indépendantes T et V 


Si nous considérons que 7 et V sont les variables indépendantes, alors : 


OU 


au 
dU(T,V) = Gr), dT + Gr). dv 


injectons cette relation dans l’expression du 1“ principe de la thermodynamique 


ÔQ = aU — dW = dU + PdV (6.5) 
pour avoir : 
5Q = (=) dT + (=) + pl dV 6.6 
e = (JT 7 OV) + 6) 
Cy l 
avec 


Cy = Gr), = eo) : capacité thermique (ou calorifique) à volume constant, [kJ/K] 


_ (au _ (82) . saat 
[= ( ov) ae = Gr), : chaleur latente de dilatation, [kJ/m' | 


? C’est des coefficients qui permettent la mesure de la chaleur mise en jeu par un système 
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> Variables indépendantes T et P 


dH(T, P) = (=) dT + (=) dP 6.7 
77 NƏT)p OP) + sae, 
Comme par définition H = U + PV , alors 

dH = dU + d(PV) = dU + PdV +VdP = 6Q+VdP (6.8) 


l’égalité entre (6.7) et (6.8) permet d’écrire : 


so = (Gr), «7 + (Gr), -V | 


Nm‘ 


Cp h 


OH ô ° 2 G \ : 

Cp = (=) = (©) : capacité thermique a pression constante, [kJ/K] 
OT/ p aT / p 

ðH 


_ (9H\ _y _ f#Q\ . 
— ( a V = G), : chaleur latente de compression, [kJ/kPa] 


> Variables indépendantes P et V 
Utilisons les définitions des différentielles totales dU et dH ainsi que les relations (6.5) et 


(6.8) : 


dU = (5) dP + (5 dV et ai= (=) dP + (=) dv 
~ \aP/y AV) p -Py OV jp 
5Q = dU + PdV 50 = dH — VdP 
y y 
50 = (5) dP + (E) + P| dV 50 = (=) v| dP + (=) dv 
v= 0P}y AV) p aa ðP )y AV) p 
—— = rad 


avec 
5Q | 
A= (©) : chaleur isochore de compression, [kJ/kPa] 
dP/y 


u = (2) : chaleur isobare de dilatation, [kJ/m°] 
aV/ p 


Remarque: Les deux coefficients calorimétriques À et u sont peu utilisés en 


thermodynamique, contrairement aux capacités calorifiques Cy et Cp qui sont très connues. 


6.3 La chaleur spécifique (à P ou à V constant) 
C’est la chaleur nécessaire pour faire varier la température, d’une unité de la matière 


(solide liquide ou gaz), de un degré (à P ou V constant, Fig.6.2). 
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m= 1 kg 
AT =1°C=1K Figure 6.2: Schéma montrant la 


définition de la chaleur spécifique 


= 5 kJ/kg °C 





On compte généralement deux types de chaleurs spécifiques: 

e La chaleur spécifique massique notée c (rapportée à la masse m d’une substance), en 
[kJ/kg.K.] = [kJ/kg.°C] 

e La chaleur spécifique molaire notée Cc (rapportée au nombre de molen d’une 


substance), en [kJ/kmol.K] = [kJ/kmol.°C]. 


Ces chaleurs spécifiques sont liées à la capacité thermique (définie précédemment, à pression 


ou à volume constant) de la façon suivante : 
C=mc=nc > c= 
v apression constante : 


Cy p = N. Cp => Ch = y 


Cp, Cp : respectivement, la chaleur spécifique massique et molaire a pression constante 


= m.c 


. OH ðh _ an 
comme par définition : Cp = |—1) alors cp =(—) etc» =(— 
P P OT P P OT P P P 


h,h : respectivement, l’enthalpie massique et molaire 


v avolume constant : 
s Cy 
Cy = M. Cy = N. Cy => Cy = 7 
M 
Cy, Cy : respectivement la chaleur spécifique massique et molaire à volume constant 
OT aT 


PPT OU ð B a 
comme par definition : Cy = (=) alors cy = (=) et Cy = ( =) 
V V V 


u,u : respectivement, l’énergie interne massique et molaire 
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6.3.1 Chaleurs spécifiques de gaz parfaits 

La chaleur spécifique est une propriété intensive de la matière qui varie essentiellement 
avec la température et l’atomicité du gaz (Fig.6.3)". La connaissance de cette variation (de 
nature non linéaire) est nécessaire pour effectuer certains calculs thermodynamiques tels que : 
la détermination de la variation de l’énergie interne (enthalpie, entropie) ou de la quantité de 
chaleur mise en jeu par un système thermodynamique évoluant entre des températures 
différentes. Mais selon les moyens de calculs (manuels ou informatiques) et le degré de 
précision souhaité, nous pouvons considérer des chaleurs spécifiques constantes ou variables. 


C 7 
kJ/kmol.K 


Figure 6.3 : Chaleur spécifique 


molaire à pression constantes 
Ar, He, Ne, Kr, Xe, Rn 


de quelques gaz parfaits 





1000 2000 3000 
Temperature, K 


> Chaleurs spécifiques constantes (valeurs approchées) 

Pour des calculs thermotechniques, qui n’exigent pas une grande précision, on peut 
considérer que la chaleur spécifique est seulement fonction de l’atomocité des gaz (donc 
indépendante de la température). Ces valeurs sont données par le tableau 6.1 (ce qui est 


valable pour AT < 100°C) (R= 8.314 kJ/kmol.K) 


Tableau 6.1 : Chaleurs spécifiques des gaz parfaits en fonction de l’atomicité du gaz 
Chaleur spécifique molaire | Chaleur spécifique massique 


Atomicité du gaz Cp 
[kJ/kmol.K] [kJ/kmol.K] nie K] Lie K] 


Monoatomique 





4 Di : P . Z ey à 
Les chaleurs spécifiques des gaz parfaits sont des fonctions croissantes de la température et de l’atomicité des 
gaz. Mais pour une substance quelconque ces propriétés peuvent varier avec aussi avec la pression ou le volume 

spécifique (inverse de la masse volumique). 
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Diatomique 
(H>, O», alr.. ) 


polyatomique 
(H20, COs,...) 





Pour des calculs un peu plus précis et lorsque les écarts de température sont de l’ordre de 
quelques centaines de degré, on peut considérer que les chaleurs spécifiques varient 


linéairement avec la température sur un écart AT = T, —T;. Dans ce cas, on utilise des 
ee ae x Tə | 
chaleurs spécifiques moyennes estimées sur AT, notées Cmlr, (ces chaleurs sont massiques ou 


molaires, à P ou V constants) . On peut les définir des façons suivantes : 
e la moyenne arithmétique des chaleurs spécifiques a T, et T; 


Tə _ Cp(Tz) + Cp(74) T, _ Cy (Tz) + cy (74) 
Cp, D Cymr a (6.9a) 


2 
Cp(T2) + p(T) T  Cy(T2) + Cy(T,) 


7 D = E 
pm |r = D ) ymy = D (6.9b) 


les valeurs de Cp, Cy, Cp ou Cy à une température T (T; ou T2) sont lues à partir du 


tableau 6.3 


e la chaleur spécifique à la température moyenne Tm = (T> + T,)/2 
T. T. 
Cpm |y = Com dm) ) Cv m| y aa CymUm) (6.10a) 
D: S _ Tz L 
pmr, = Ch) J Cv ml, = CymUm) (6.10b) 
e ou encore 


Tz Cp m(T2). Log Cp m(1). Ti Tz Cy m(T2). T; — Cy m(Ti). Ti 


Pml y, ne À vml, rer (6.11a) 
= ie = Cp m(T2). T2 = Cp m(T1). Ti 5 ie = Cy,m(12).T2 = Cy m(T1). Ti (6 11b) 
P,m Ti T» = Ti J V,m T T; = A _ = . 


Ces deux approximations ne sont pas utilisées pour le calcul de AS 


Notons que les Cpm(T) et Cym(T) (ou Cpm(T) et Cym(T)) sont données dans des tableaux 


sous forme de valeurs directes ou sous forme de polynôme de premier degré (Tableau 6.2). 


Donc si par exemple on veut déterminer la variation Au ou Ah entre deux températures 
Ti etT2 , on écrit : 
T2 


T: T: 
Au = | cy(T) dT = CV m . .AT = Cv mly (T> = Ti) 
T1 
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T2 
T- T: 
Ah = f Cp(T) dT = cpm x ATS Cpm ly (D — T;) 


T1 


Tz Tz 2 X\ ° 9 . 7 2 
avec, C et C calculées à partir de l’une des relations précédentes 
V,m T P,m IA 


> Chaleurs spécifiques variables (valeurs exactes) 
La variation réelle (exacte) de la chaleur spécifique avec la température est disponible 
dans la littérature sous forme de graphe (Fig. 6.3), de valeurs tabulées (tableau 6.3) ou de 


polynômes (du 3° degré en T et plus) tabulées (tableau6.4). 


6.3.2 Chaleurs spécifiques des liquides et solides 
Pour un liquide ou un solide, la chaleur spécifique à pression constante est égale à la 
chaleur spécifique à volume constant: 
Cp = Cy =C 
ces propriétés thermodynamique ne dépendent que de la température (augmente avec Ty 


comme dans le cas des gaz parfaits. 


6.3.3 Relations entre les chaleurs spécifiques d’un gaz parfait 
Dans le cas de gaz parfais, l’ énergie interne massique u et l’enthalpie massique h ne 
dépendent que de T, ce qui nous permet d’écrire : 


= (=) _ du : = (=) _ dh = 
Vv \ar), ar © PT dT 


d(h—u)  d(Pv) 
dT dT 





Cp — Cy = 


comme le gaz est parfait (Pv = rT), alors : 


Cp— Cy =r (6.12) 


En multipliant par la masse molaire M du gaz, nous obtenons : 


Soit y la constante adiabatique d’un gaz, définie par : 


Cp Cp 
y = — | ou — = Cp =Y.Cy 
Cy Cy 


| Excepté pour l’eau dont la chaleur spécifique diminue jusqu’à 35°C puis augmente 
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En injectant cette relation dans (6.12) on obtient : 
r : V.T 
Cy=— et Cp =—— 
” -1 © (y-1) 
En multipliant ces relations par M, nous aurons : 


R y.R 


v= GD &= Ga 


Exemple : de lair a 300K et 200kPa est chauffé a pression constante jusqu’à une 
température de 600K. Calculer la variation de l'énergie interne massique en 
utilisant différentes chaleurs spécifiques : (a) constante, (b) moyenne et (c) variable. Quelle 
erreur de calcul va-t-on commettre dans chacun des cas si on connait les valeurs exactes (de 
mesure) de u : u (300K) = 214.07 kJ/kg et u (600K) = 434.78 kJ/kg (Mair = 28.97 kg/kmol) 


Solution : 
AUpyq = 437.78 — 214.07 = 220.71 kJ/kg 


(a) On lit la valeur de la chaleur spécifique massique constante a partir du tableau 6.1 


cy = 21 = 0.717 kf /kg.K . Donc : 


600K 
Auta) = f Cy dT = cy AT = 215.10 kJ/kg 
300K 
Boe AN 
ETT(a) = mii A te 


AUexa 
(b) on utilise le tableau 6.2 pour calculer la valeur moyenne de la chaleur spécifique à la 


température moyenne donnée par la relation (6.10a). : 


600K 


600K 
AU(p) = | Cy dT = T (T- = T4) = 217.5 kJ/kg 
300K 
avec :Cy m a = Cy m(Tm = 450K) 


= 0.7084 + 0.00009349 (t = 177°C) = 0.725 kJ /kg.°C 
Rappelons que la température t, dans les polynômes donnés par le tableau6.2, doit être en 
degré Celsius °C et non en Kelvin. 


AUexa —AU(p) 


ETT(p) = = 1.45% 


AUexa 
(c) le tableau 6.4 donne l’expression de la chaleur spécifique molaire à pression constante 


(Cp). Pour avoir la chaleur spécifique massique à volume constant (cy) on écrit : 
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Cp(T) 
M — r 


R 
= g (a+ BT +yT* + 6T* +eT*) -r 


Cy(T) = CT) -r = 





= r((a — 1) + PT +yT? + ÔT* + ET*) 


donc : 
600K 600K 
Auco) = f cy(T)dT =r | ((a-—1)+ BT +yT* + 6T? + eT*)) aT 
300K 300K 
B y 5 600K 
= SO ee ee 
r((a ) +5 +3 Lab +e nook 


pour l’air: a = 3.653 , B = —1.337 x 107? , y = 3.294 x 1078 ,6 = —-1.913 x 10 ° 


e = 0.2763 x 1071? 


Auco = 220.69 kJ/kg 


AUexa —AU(c) 


eTT(c) = = 0.009% 


Alexa 


Tableau 6.2: Variation linéaire des Chaleurs spécifiques massiques moyennes a P ou V 
constants des gaz parfaits 


Cpm(t) Cy,m(t) 


0.9952 + 0.00009349 t 0.7084 + 0.00009349 t 
H2O 1.84012 + 0.00029278 t 1.3783 + 0.00029278 t 


Intervalle de température 0 à 1500°C 


14.2494 + 0.00059574 t 10.1241 + 0.00059574 t 
0.8725 + 0.00024053 t 0.6837 + 0.00024053 t 
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Tableau 6.3 : Chaleurs spécifiques massiques des gaz parfaits [kJ/kg.K] 


ete tete tet eto La Le 


| 1.030 | 0.743 | 1.40 ns | oess | 1395 | 300 
e | 1.041 | 0.744 | 1 28 | 189 | 350 


1.110 36: 1.043 783 431 750 
1.121 l. 054 À g $00 
Temp. 


s 523 288 ses ggg 7] 


Ae 


= 282)" 


i 
= 


BS £22 28% Lt 


| 1.163 | 0.866 14.822 10.6 1.385 ox) 
1.185 | 14.983 10.85 | 380 1000 


Source: Adapted from K. Wark, Thermodynamics, 4h ed., McGraw-Hill. New York, 1983, as based on “Tables 
of Thermal Properties of Gases,” NBS Circular 564, 1955. 
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Tableau 6.4 : Variation de la chaleur spécifique molaire à P constante des gaz parfaits 


Cp = R(at+ BT + yT? + 8T? +ET*) 
T is in K, equations valid from 300 to 1000 K 


Gas a B x 10° y x 10° & x 10° e x 10° 
CO 3.710 — 1.619 4.692 —2.032 0.240 

H; 3.057 2677 —3.810 5.521 —1.812 

H,0 4070 — 1.108 4.152 — 19064 0.807 

0; 3.626 — 1.878 7.055 — 6.764 2.156 

N; 3.675 — 1.208 2.324 — 0.632 0.226 

Air 3.653 = 1.337 4.294 — 1.913 0.2763 
50, 3.267 5.324 0.684 — 3.281 2.559 

CH, 3.826 —3.979 24.558 22733 6.963 

CH: 1.410 19.057 — 24.501 16.391 —4,135 

CHa 1.426 11.383 7.989 =16.254 6.749 

Monatomic 

gases“ 25 0) 0 0 0 


“For monatomic gases, such as He, Ne, and Ar, 7, is constant over a wide temperature range and is very nearly 


equal to 5/2 R 


Source: Adapted from K. Wark, Thermodynamucs, 4th ed. McGraw-Hill, New York, 1983, as based on NASA 
SP-273, US. Government Punting Office, Washington, DC, 1971. 
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Chapitre 7 : 


Le deuxième principe de la thermodynamique 


7.1 Insuffisance de 1°” principe de la thermodynamique 

Soit un corps chaud à la température T; placé dans un milieu ambiant chaud Tọ (Fig.7.1). 
Après un certain temps, le corps va spontanément s’échauffé (cas (a)). Cette transformation 
satisfait parfaitement le premier principe de la thermodynamique (ou principe de conservation 
de l’énergie) puisque la quantité de chaleur Q cédée par le milieu ambiant chaud est égale à la 
chaleur gagnée par le corps froid. Imaginons, à présent, le refroidissement supplémentaire du 
corps froid a partir du milieu ambiant chaud (cas (b)). Même si cette transformation est 
impossible elle respecte toujours le premier principe de la thermodynamique : la chaleur reçue 
par le milieu chaud est égale a la chaleur cédée par le corps froid. Donc, il est clair que 
l’évolution d’une transformation est possible dans une direction et impossible dans l’autre. 
Mais, le premier principe est incapable de nous fournir cette information. Pour compléter 
cette insuffisance, un second principe de la thermodynamique a été établi. Ce dernier fait 
intervenir une nouvelle fonction d’état appelée entropie et notée S. 


20°C 





5°C 


O T \ a! 


(a) Echauffement (b) Refroidissement 











Figure 7.1 : Exemple de transformations thermodynamiques 
(a) irréversible, (b) impossible 


7.2 Notions de bilan d’entropie 
Pour mieux comprendre comment s’effectue un bilan d’entropie, nous allons tout d’abord 

envisagé la situation suivante : soit un petit village d’un certain nombre d’habitants N. Ce 
nombre d’habitants peut varier pour deux raisons possibles : 

- l’arrivée ou le départ de certaines personnes étrangères ; 

- la naissance ou le décès de personnes à l’intérieur du village. 
Donc, pour estimer à tout moment la variation de la population du village, on effectue le bilan 
suivant : 

AN = AN, + AN; 
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avec : 
AN : variation totale de la population; 


AN, : variation externe de la population, due aux départs et arrivées 
AWN; : variation interne de la population, due aux naissances et décès 
Remarque : 
Lorsque AW, = 0 (aucun départ ni arrivée) et AN; = 0 (aucune naissance ni décès) alors 


AN = 0. Cette situation correspond à un état de conservation de la population du village. 


On peut tout à fait effectuer un bilan analogue à un système thermodynamique qui échange 
de l’énergie thermique (chaleur) avec le milieu extérieur. Dans ce cas c’est l’entropie qui va 


caractériser la variation de l’énergie du système 


7.2.1 Bilan pour un système fermé 


Pour un système thermodynamique fermé subissant une transformation élémentaire, on 


écrit : 
ÔQe 
dS = Sse T Ssi = —— + Ssi (7.1) 
=x ‘i 
ÔQe 
T 
avec : 
dS : variation élémentaire de l’entropie du système; 
Sse: entropie transférée due aux échanges de chaleur (6Q, )! avec le milieu extérieur; 
So; : entropie créée à l’intérieur du système (ou terme source) due aux irréversibilités 
internes. 
Remarques 


e Sse > 0 (échauffement du système), < 0 (refroidissement), = O (adiabatique) 

e Sco; > O (transformation irréversible), = O (réversible). Le terme source de l’entropie 
n’est JAMAIS négatif 

e Pour un système adiabatique (Sse = 0) dS = Ssi (= 0), ceci correspond à l'un des 
énoncés du deuxième principe de la thermodynamique qui stipule que : l’entropie d’un 


système adiabatique ne diminue jamais 


En intégrant la relation (7.1) le long d’une transformation réversible et irréversible 


évoluant entre les mêmes états d’équilibres initial et final (Fig.7.2), on peut écrire : 


* 6Q, > 0 (échauffement du système), 6Q, < 0 (refroidissement) et 6Q, = 0 (système adiabatique) 
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(2) (2) 
ASe, = Sy — S, = | dS = | — (7.2) 
(1) (1) 
(2) (2) 
AS; = S, — Sy = fas=]( Cert + Ssi) (7.3) 
(1) (1) 
Comme l’entropie est une fonction d’ état, alors: 
(2) (2) 
BSirr = ASrev = | “ere > | er (7.4) 
(1) (1) 


AS=S,-5,=0.4 kI/K 






Transformation 
irreversible 


\ 2 
T 
oe 
F 


b 







Figure 7.2 : variation de l’entropie 
pour une transformation réversible 
et irréversible. 


Transformation! 


| 
| 
| 
| 
| 
| 
| 
| 
l 
| reversible | 


0.3 0.7 5, KIK 


7.2.2 Bilan pour un système ouvert 
Dans le cas d’un système ouvert (Fig.7.3) il faut en plus tenir compte de l’entropie 


transférée par la masse qui transite (rentre et sort) du système. On écrit alors : 


as = y& E tms- Dm s + [ss [kJ/K] (75) 
v 
Variation 
d'entropiedu T. , — terme source 
entropie 
système ouvert ae ape U PENON 
transferée transferée total 
par échange par échange 
de chaleur de matière 
| . 
| = 
| Sj kTkgK] 
l 
Système ouvert | 


N Figure 7.3 : Bilan d’entropie d’un 
4 QO. système ouvert 





en se rapportant au temps, le bilan s’écrit : 
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Q Ssi 
ii X m, Si — Xn Si t T [kW/K] (7.6) 
i j 
avec, Q, = 6Q,/dt et m = dm/dt 


> Régime permanent (dS/dt=0) : 


DEL “Dist | ae (7.7) 


> Régime permanent, un seul flux de matière : 
oo 
0 = > Cr + MS, — S2) + — (7.8) 
= Tk dt 


> Régime permanent, un seul flux de matière, adiabatique (Fig.7.4): 


Ssi 


dt 
> Régime permanent, un seul flux de matière, adiabatique, réversible (Fig.7.4): 


0 _ m(s: = S>) => S1 = S2 (7.10) 


0 = m(s,; — S2) + (7.9) 


Éi Pi; Si 
| | _— | — | 


4 Pat i i 

, l me —— 

| irreversible [IN | reversible ara 

| 4 4 
Le) 





isolation = lp m. \ 4 
gi | H 
thermique 


Figure 7.4 : Schéma montrant la variation de l’ entropie pour 
une détente adiabatique réversible et irréversible 


7.3 Variation de l’entropie d’un gaz parfait 
Nous allons, dans ce paragraphe, calculer la variation de l’entropie AS d’un gaz parfait 


subissant une transformation réversible dans un système fermé. 


e AS en fonction de T et V : 


premier principe dU = —PdV +Tds 


SW 5Q = dT dV 71] 

gaz parfait dU = m cy dT = $= m| cy +r (7.11) 
i re 
y 


en intégrant cette relation entre les états (1) et (2) et en supposant Cy constante : 
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AS = | in + rin 
=m|cv lny trny 


Tə 
= n| yn 


A+R ni [kJ /K] 


e AS en fonction de TZ et P : 


premier principe dH =VdP+Tds 


dS = | cp z = (7.12) 
m| c — r — : 
gaz parfait dH = m cp dT = eye P 
y= mrT 
HORS 

en intégrant cette relation entre les états (1) et (2) et en supposant Cp constante : 

T> P, 

AS = m | cp morin 
P, 


T 
= n| Zp nz- R n2 [kJ /K] 
1 


1 
e AS en fonction de Pet V: 


PV=mrT = InP +]InV = InT +In(mr) > 


dP dV dT 
P V T 
Substituant cette relation dans (7.12) pour avoir : 
dP dV dV dP 
dS = m| cp (= +£) -rF]= m| cp + Cy — (7.13) 
V P 
en intégrant : 
V» P> 
AS = m | cri + Cy In 
P, 


Be E 
= n| Gin + Cy m 
1 1 


7.3.1 Cas particuliers de transformations réversibles 








e Transformation isobare : dS = ôQ dH _ mcpaT de wee 
oe y aA 
TARE ôQ dU  mcyaT 
e Transformation isochore : dS = — = — = NS = l 
g T meyln g- 
oe ôQ —óôW  PdV dV V> 
e Transformation isotherme : ds = — = — = — = mr — => NS = mrin < 
T T T 1 
| | | dQ 
e Transformation adiabatique : dS = 7 = 0 = S = cste 
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e Transformation polytropique : 


TV* t=cste => InT + (k — 1)InV = cste => 
à n =o ne dT 1 
T V V2 T (k—1) 


Injectons cette égalité dans la relation (7.11) pour avoir : 


E r dT cy(y—1)\aT aT 
aS =m (ey Gay) = (yao) eee 


7.3.2 Diagramme entropique (T-S) 
Dans un diagramme entropique, la surface en dessous d’une courbe correspond a la 


quantité de chaleur mise en jeu le long d’une transformation thermodynamique (Fig.7.5) 








Transformation 
reversible 


aa = T'as = ôg 
/ 


= | 
Ares =| TdS =Q! Figure 7.5 : Diagramme entropique 





Dance ce diagramme une isotherme est représentée par une droite horizontale et une 
isentropique par une verticale, mais pour une isobare, isochore ou une polytropique il faut 
chercher les équations des courbes qui caractérisent ces transformations. 


e allure d’une transformation isobare : 


T S 
ac aT aT dS | T S=- _, 
= — = — Z — — = — 
ee T T MCp "T MCp 

16 So 
(0) 
T =7T,.e "CP/ = branche exponentielle 





l’isobare est une branche exponentielle croissante de pente 
aT T ei 
dS mcp 
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e allure d’une transformation isochore : 


dT (2—0) 
dS = mcy F => T =7T,.e Mcv/ = branche exponentielle 





l’isochore est une branche exponentielle croissante de pente : 


dT 
dS mcy 
pente isochore Cp 


=—=y>r>0 => 
pente isobare Cy £ 


>0 


dans un diagramme entropique une isochore est plus raide qu'une isobare 


e l’allure d’une transformation polytropique : 





dans ce cas la pente de la courbe peut être positif ou négatif : 


dT T k—y 
Go we Car VE] 





ko+c0 T 





Figure 7.6 : Représentation de transformations thermodynamiques 


dans un diagramme (P-V) et (T-s) 
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7.4 Enoncés du deuxième principe de la thermodynamique 
74.1 Définitions 
e Source de Chaleur: Une source de chaleur à la température T est un système 
thermodynamique suffisamment grand pour que sa température reste constante 
lorsqu’il cède ou reçoit de la chaleur, exemple : l’atmosphère, un fleuve, un lac, un 


four réglé à une température...etc 


e Cycle monotherme : cycle d’un système qui échange de la chaleur avec une seule 


source de chaleur. 


e Cycle ditherme : cycle d’un système qui échange de la chaleur avec deux sources de 


chaleur (chaude et froide). 


74.2 Enoncé de Kelvin-Planck 

Il n’existe pas de machine qui soustrait de la chaleur à une seule source de chaleur pour 
fournir une quantité équivalente de travail (Fig. 7.7). Autrement dit une machine ne peut 
fonctionner selon un cycle monotherme, ou encore qu’il est impossible qu’une machine 


possède un rendement (gain d’énergie/ dépense d’énergie) égal à un. 


Ca 
= 


= wre 


Source de chaleur . . . 
Impossible sauf si la machine 






| échange de la chaleur avec une 
Q = 80 kW deuxième source de chaleur 


W = 80 kW 





Machine 


Figure 7.7 : Schéma d’une machine qui ne respecte 


pas l’énoncé de Kelvin- Planck 


74.3 Enoncé de Clausius 
Il n’existe pas de machine qui transforme spontanément de la chaleur d’une source froide 


à une source chaude (Fig.7.8) 
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Source chaude 






Oc=5kW 
| Impossible sauf si on fournit à la 


N machine du travail 
Machine 


Sourcetroide 


Figure 7.8 : Schéma d’une machine qui ne respecte 
pas l’énoncé de Clausius 


7.4.4 Enoncé de Carnot 
Pour faire fonctionner une machine, celle-c1 doit échanger du travail, avec le milieu 


extérieur, et de la chaleur avec au moins deux sources de chaleur (chaude et froide) 


Source chaude _ 


Qc = 80 kW 





Possible 


Source froide 


Figure 7.8 : Schéma d’une machine qui respecte l’énoncé de Carnot 
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Chapitre 8 : 


Cycles des moteurs thermiques 


8.1 Définitions des moteurs thermiques 

Un moteur thermique est une machine qui transforme de l’énergie thermique en une énergie 
mécanique : 1l reçoit, en un cycle (par l’intermédiaire d’un fluide moteur), une quantité de 
chaleur Qc de la source chaude (a la température Tc) et rejette une quantité de chaleur Q;a la 
source froide (à la température Ty) (Fig. 8.1). Une partie de la quantité de chaleur Qc reçue est 


transformée en un travail W que le moteur fourni au milieu extérieur. 


Qc>=0 


SYSTEME 


FLUIDE 
THERMIQUE 


| | TRAVAIL 
r tL W<Ù 
a = 


Figure 8.1 Schéma de fonctionnement d’un moteur thermique 





8.2 Le rendement thermique d’un cycle moteur 

Le rendement thermique d’un cycle est le rapport entre l’énergie utile (ou valorisable) et 
l’énergie coûteuse qu’il a fallu fournir. Autrement dit, le rapport de ce qu’on veut (W) sur ce 
qu’il faut dépenser pour l’avoir (Qc). 


gain énergie utile  [Woyce 





Nth <1 





dépense énergie coûteuse Qe 


Remarques 
= Les moteurs thermiques sont des machines motrices (qui produisent du travail), leurs cycles 
a i 
ne sont donc parcourus que dans le sens des aiguilles d’une montre. 
= Dans tout ce qui va suivre, les fluides qui vont décrire les cycles sont supposés être des gaz 


parfaits. 
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8.3 Le cycle de Carnot 

On appelle cycle de Carnot, du nom de l’ingénieur français Carnot, le cycle idéal constitué 
de quatre transformations réversibles dont deux isothermes (de températures égales a celles des 
sources) et deux adiabatiques (isentropiques). La représentation de ce cycle dans le diagramme 


de Clapeyron (P,v) ou le diagramme entropique (T,s) est illustrée par la Figure8.2. 





Figure 8.2: Le cycle de Carnot 


a) Les différentes étapes du cycle 
Compression 1sentropique (Q12 =0) 
Wis = W;, = Au, = C,(Te — Tr) 
Détente isotherme (Au,3 = 0) 
W3 = —Qc = —Tc(S3 — S2) avec (S3= S4 et So= Sı) 
Détente isentropique (Q34=0) 
Wis = W3, = Aus, = Cy (Ty — Te) 
Compression isotherme (Au,, = 0) 


Wii = —Q¢ = —T; (Sy — Sy) 





Le bilan du cycle 


Wutile = cycle — Wiz + Wo3 + W34 + Way = — (Qc + Qr) 
Energie utile 
= (Si — S4)(Te — Tf) 


Energie perdue Qr = T (Sy — S4) 
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Rendement 


thermique du cycle 





Théorème 1 de Carnot : 
Le rendement thermique d’un moteur fonctionnant selon le cycle de Carnot est indépendant 


de la nature du fluide moteur et ne dépend que des températures des sources chaude et froide. 


Théorème 2 de Carnot : 
Le rendement d’un moteur réel est toujours inférieur au rendement de Carnot qui constitue 
la limite théorique maximale impossible à dépasser (n1 à atteindre) quel que soit la perfection 


technique de la machine. 


Remarque : 

Un moteur qui fonctionne selon le cycle de Carnot a un rendement maximal mais présente 
l’inconvénient d’être très compliqué (donc très coûteux) à réaliser, pour des raisons techniques 
et non théoriques (isothermes difficiles à assurer, moteur encombrant,...etc). Pour cela, on 
préfère fabriquer des moteurs ayant un rendement plus faible mais plus simples (moins chers) à 


mettre en œuvre. 


8.4 Le cycle d’Otto ou de Beau de rochas (1862) 

C’est le cycle théorique selon lequel fonctionne un moteur à essence à 4 temps. L’apport de 
chaleur (Qc) se fait à volume constant (et non pas à T constante). Ce cycle se compose de deux 
transformations isochores (2-3 et 4-1) et de deux transformations adiabatiques (1-2 et 3-4) que 
l’on supposera réversibles. L’allure du cycle d’Otto dans le diagramme de Clapeyron (P,v) ou 


le diagramme entropique (T,s) est illustrée sur la Figure 8.3. 





Figure 8.3 : Le cycle d’Otto 
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a) Les différentes étapes du cycle 
Admission : aspiration du mélange air + essence (ce n’est pas une transformation au 
sens thermodynamique) 
Compression isentropique (Q72 =0) 


Wi; = Wiz = Aus = C,(T2 — T,)> 0 


T, = (a tT (équation caractéristique d’une transformationisentropique) 


Apport de chaleur (combustion) à volume constant (W:3= 0) 

Qc = Aug3 = Cols — T2) = CyT2(T3/T2 — 1)> 0 

Détente isentropique (Q34=0) 

Wis = Was = Augy = Co(Ts — T3)< 0 

T; = (a’1)T, (équation caractéristique d’une transformation isentropique) 
Refroidissement isochore (W4; = 0) 


Qs = Aug, = Cols — Ty) = —C,7T,(T,/T, — 1)< 0 
1-0 Echappement : évacuation des gaz brulés 


y : la constante adiabatique du fluide moteur (= C,/C,) 





a : le taux de compression (= V,/V; =V4/V3> 1), c’est un paramètre de construction du cylindre 


moteur 
Remarque : 


= Les transformations citées ci-dessus vont avoir lieu à l’intérieur du cylindre moteur, elles 


correspondent aux 4 temps suivants : 


aS 























e ee 


| Temps 2 Temps 3 Temps 4 Temps 
(Admission) (Compression) (Combustion — Détente) (Refroidissement — Echappement) 
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Un seul cylindre moteur ne peut pas produire un mouvement de rotation (du vilebrequin) 


continue, ce qui nécessite la présence de trois autres cylindres qui travaillent chacun sur un 


temps différent : 





Vilberequin 
Le bilan du cycle 
Wutite = Weycte = Wiz + W34 = —(Qc + QF) 
s 
Rendement To) = 1 _ +i -jai 





thermique du cycle 





Remarque : 
Le rendement thermique du cycle d’Otto ne dépend que du taux de compression (a)et de la 


constante adiabatique (y).Pour un fluide moteur donné, plus (a) est grand, plus le rendement 


thermique est élevé (Figure 8.4) : 


OS 








C,H, 
Typical 
compression 
¿Æ ratios for Figure 8.4 : Le rendement thermique du 
gasoline cycle Otto en fonction du taux de 


ee ete compression a 
10 12 


là 

+ 
-a 
Pr 
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8.5 Le cycle Diésel 

C’est le cycle théorique selon lequel fonctionne un moteur Diesel à 4 temps. L’apport de 
chaleur (Qc) se fait à pression constante. Ce cycle se compose d’une transformation isobare (2- 
3), une transformation isochore (4-1) et de deux transformations adiabatiques (1-2 et 3-4) que 
l’on supposera réversibles. L’allure de ce cycle dans le diagramme de Clapeyron (P,v) ou le 


diagramme entropique (T, s) est illustrée sur la Figure 8.5 





Figure 8.5 : Représentation du cycle Diesel sur le diagramme (P, V) et (T,S) 
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a) Les différentes étapes du cycle 


0-1 Admission d’air seul 


Compression isentropique (Q72 =0) 


Wi; = Wy2 = Aus = C,(T2 — T,)> 0 


T> = (a "NT, (équation caractéristique d’une transformation isentropique) 


Apport de chaleur (combustion) à pression constante: injection de gasoil 

Qc = Ahg3 = C,(T3 — T2) = CyT2(T3/T2 — 1)> 0 

W23 = Auz3 — Qc= (C,-C, (T3 — T2)< 0 

T3/Tz = V3/V2 = b (équation caractéristique d’une transformation isobare) 
Détente isentropique (Q34-0) 

Wis = Was = Aug, = C,(Ty — T3)< 0 

T; = (b/a)'-’T, (équation caractéristique d’une transformation isentropique) 
Refroidissement isochore (W4; = 0) 

Qf = Aug, = Cols — T4) = —-CT1(T4/T1 — 1)< 0 


1-0 Echappement : évacuation des gaz brulés 





a : le taux de compression (= Vı/V2 > 1) 


b : le taux d’ injection (= V3 /V2> 1) 


Les transformations qui composent le cycle Diésel correspondent aux 4 temps suivants : 


aT 


(4) 














(Admission) (Compression) (Combustion — Détente) (Refroidissement — Echappement) 
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Le bilan du cycle 


Energie utile Wutile = — cycle — = W2 + W34 + W3 = = —(Qc T Qr) 
e aa 


Qr 
Qc 


Rendement 


Weycte 


=1+—4 
Qc 








ang 





thermique du cycle a | 





Application numérique : P;= 1bar, T7=55°C, y =1,4 et b —2 





Les valeurs usuelles du taux de compression a sont entre 14 et 20. Celles du rapport du taux 
d’injection b entre 2 et 2.5. Le rendement thermique est usuellement compris entre 60 et 65%. 


L’évolution du rendement thermique du cycle Diésel en fonction de a et b est donnée par la 


figure 8.6 
T 
b= L 

0.6 _ 

- 0.5 à 

4 0.4 

= | Typical 
0.3 compression 
02? ratios for diesel 
o | engines 
O.1- | Figure 8.6: Rendement thermique du 


cycle Diésel en fonction du taux de 
6 8 10 12 14 1618202224 compression a et d’injection b 
Compression ratio, @ 


[= 
a 
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Remarque : 

Dans le cas du moteur à essence (Fig. 8.7a), qui fonctionne selon lecycle Otto, l’apport de 
chaleur dans le cyclindre moteur se fait à l’aide d’une éteincelle provoquée par une bougie 
d'allumage. Dans le cas du moteur Diésel (Fig. 8.7b),l’apport de chaleur se fait par l’injection 
de gasoil à travers un injecteur. 


Bougie Eteincelle 


Mélange 
Air-essence 





(a) Moteur essence (b) Moteur Diesel 


Figure 8.7 :Schéma d’un cylindre moteur 


8.6 Le cycle mixte de Sabathé 


C’est un cycle d’Otto (Beau de Rochas) combiné avec un cycle Diésel (Fig. 8.8). L’apport 


de chaleur (combustion) se fait de façon isochore puis isobare. 





Figure 8.8: Le cycle mixte de Sabathé 
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a) Les différentes étapes du cycle 


Admission. 


Compression isentropique (Q72 =0) 
Wis = Wiz = Aus = C,(T2 — T,)> 0 
= (a’~1)T, (équation caractéristique d’une transformation isentropique) 
Apport de chaleur (combustion) à volume constant (W23= 0) 
Qc1 = Aug3 = C,(T3 — T2) = C,T2(T3/T2z — 1)> 0 


T3/T> = P3/P, =r (équation caractéristique d’une transformation isobare) 


Apport de chaleur (combustion) a pression constante 


Qc2 = Ah = C,(T, — T3) = CpT2(T4/T2 — T3/T2)> 0 


W34 = Au34 — Qc2 = (Cy-Cy (Ty, — T3)< 0 

Détente isentropique (Q45=0) 

Wis = Was = Augs = C,(Ts — Ta)< 0 

T, = (V,/V=)*-’T. = (b/a) YT; (équation caractéristique d’une transformation 
isentropique ) 

Refroidissement 1sochore (Ws; = 0) 

Qs = Aus, = CT — Ts) = —C,T,(Ts/T, — 1)< 0 

T/T, = P/P, = rb’ 


1-0 Echappement : évacuation des gaz brulés 


a : le taux de compression (= Vj/V2 > 1) 





b : le taux d’injection (=V4/V3> 1) 


r : le rapport de pression (= P3/P;> 1) 


Le bilan du cycle 


Energie utile Wutile = — +Y cycle — = W2 + W34 + Was = —(Qc1 + Qc2 + Qr) 
en 


Rendement rbY —1 


a | 
r—1+yr(b-1) 


"Qc 


Weycte| _ 


Qc 





en = | 1-an] 


thermique du cycle 
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Application numérique :P7= 1bar, T;=55°C, y =1,4 et b= 2 





8.7 Le cycle de Joule-Brayton 
C’est le cycle théorique selon lequel fonctionne la plupart des turbines a gaz (Fig. 8.10). Ce 
cycle est composé de deux transformations isentropiques (1-2s) et (3-4s) et deux 


transformations isobares réversibles (2s-3) et (4s-1) (voir Fig. 8.9): 


Chambre de 
combustion 





Figure 8.9 : Schéma simplifié d’une turbine à gaz 


Qc : Puissances thermiques échangées avec la source chaude, [Kw]. 
Qr : Puissances thermiques échangées avec la source froide,[Kw]. 


Weycte : Puissance mécanique du cycle, [Kw]. 
Qc = Qc.m, Qf = Qs. met Weycte = cycle: M 


m : Débit massique du fluide moteur, [Kg/s] 
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Figure 8.10 : Le cycle réversible de Joule -Brayton 


a) Les différentes étapes du cycle 
Compression isentropique (Q72, =0) 
Wis = Wizs = Ahy2s = Cy (Tzs — T1)> 0 
To, =AY-V/NT, (équation caractéristique d’une transformation isentropique) 
Apport de chaleur (combustion) à pression constante (W23 —0, système ouvert) 
Qc = Ahgs3 = Cy(T3 — Tzs) = CpT25(T3/T25 — 1)> 0 
Détente isentropique (Q34;=0) 
Wis = Wz4s = Ah3as = Cy (Tas — T3)< 0 
Tz; =AY-VIMT, . (équation caractéristique d’une transformation isentropique) 


Refroidissement isobare (W4,; —0, système ouvert) 


Qs = Alas, = Co — Tas) = —CpT1(T4s/T1 — 1)< 0 





À : le rapport de pression (= P2/ P1) 


Le bilan du cycle 


Energie utile Wautite = a cycle — = Wizs + W345 = —(Qc + QF) 
W, 
= 1 — 


Rendement 


thermique du cycle 
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Remarque : 


Le rendement thermique du cycle augmente avec le rapport de pression A= P2/ P; (Fig. 8.11) : 


d th, Brayton 


Typical pressure 
ratios for gas- 
turbine engines 





2 10 15 20 25 
Rapport de pression, À 
Figure 8.11 : Le rendement thermique du cycle Joule Brayton 


en fonction du rapport de pression, À 


8.7.3 Le cycle de Joule-Brayton avec irréversibilités 
Supposons maintenant que le cycle de Brayton (Fig. 8.12) comporte une compression 


irréversible (adiabatique) (1-2), de rendement isentropique de compression : 


: / . C 
a travail réversible o Ws Ds — Ty pe 


lis ~ travail irréversible W° T, —T, 





et une détente irréversible (3-4), de rendement isentropique de détente : 


p travail irréversible W?  T1T,—T; 
Ni = = 7 Ci — = < 1 
5  travailréversible W?  Tis — T3 
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Qc a | f 





5 


Figure 8.12: Le cycle irréversible de Joule Brayton 


a) Les différentes étapes du cycle 
Compression adiabatique (Q2 =0) 
WC = Wiz = Ah = C,(T, — T,)> 0 
We = Wis Mis 
Apport de chaleur (combustion) à pression constante (W23 =0, système ouvert) 
Qc = Ahz3 = Cy (T3 — T2) 
Détente adiabatique (Q34=0) 
W? = Way = Ahgy = Cola — T3)< 0 
WP =nisWs 
Refroidissement isobare (W4; =0, système ouvert) 


Qs = Ahy, = Ci — T,)< 0 





b) Le bilan du cycle 


Energie utile Wutite = cycle — Wiz + W34 = —(Qc + Qr) 


Rendement Weycie| _ i+ Qf 
Qc Qc 





en = | 


thermique du cycle 
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8.7.4 Amélioration du cycle de Joule-Brayton 
Il existe plusieurs façons d’améliorer le cycle de Joule (c’est à dire augmentation du 


rendement thermique) 


8.7.4.1 Cycle à régénération 
La régénération (Fig. 8.13) consiste à récupérer la chaleur des gaz chauds à la sortie de la 
turbine, par l’intermédiaire d’un régénérateur (échangeur de chaleur), pour chauffer les gaz 


après la compression et avant la combustion (ce qui n’est possible que si Ty, >7>,) : 


=, Chambre de 
| combustion 





Figure 8.13 : Le cycle de Joule Brayton avec régénération 


a) Les différentes étapes du cycle 
Compression isentropique (Q72; —0) 
Wis = Wizs = Ahizs = Colis — T1) > 0 
Préchauffement (isobare) de I’ air dans le régénérateur (récupération de la chaleur) 


Qrég = Ahzss = Cy (Ts — Tp) > 0 


Apport de chaleur (combustion) à pression constante (W53 —0, système ouvert) 


Qc = Ahs3 = C,(T3 — Ts) > 0 


Détente isentropique (Q34,;=0) 
W = Was = Ahzas = Cp(Tas — T3)< 0 


L 
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Prérefroidissement (isobare) des gaz brulés dans le régénérateur 


Qrég = Alase = Cp(Te — Ta) < 0 


Refroidissement isobare (W6; =0, système ouvert) 


Qs = Ahe, = Cols — Te) < 0 





À : le rapport de pression (= P2/ P1) 


Le bilan du cycle 


Energie utile Wutile = — cycle — = W425 + Was = —(Qc + Qr) 
Teeme o 


Rendement Weycle 


Nth Qc 








thermique du cycle 





8.7.4.2 Cycle avec fractionnement de la détente 


Le principe est d’effectuer une surchauffe intermédiaire entre plusieurs détentes adiabatiques 


de taux de détente moindre. 























Compresseur 


T 


mi þa q e | = 
= | 





Figure 8.14 : Le cycle de Joule Brayton avec fractionnement de la détente 
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a) Les différentes étapes du cycle 
Compression isentropique (Q72; =0) 
Wis = Wizs = Ahiss = Cols — T4) > 0 
Apport de chaleur (combustion) à pression constante (W,3 = 0, système ouvert) 


Qc1 = Ahgs3 = Cy(T3 — T25) > 0 


Détente isentropique dans le 1“ étage de la turbine (Q3as = 0) 


Wis = Was = Ahzas = Cy (Tas — T3) <0 


Détente isentropique dans le 2° 


Wis = W3as = Ahzas = Cy (Tas — T3) <0 


étage de la turbine (Qasbs = 0) 


Prérefroidissement (isobare) des gaz brulés dans le régénérateur 
Qrég = Muse = Cp(Té — Tas) < 0 
Refroidissement isobare (Ws; =0, système ouvert) 


Qs = Ahey = Cy) — Te) <0 





À : le rapport de pression (= P2/ P;) 


c) Le bilan du cycle 


Energie utile Wutile = cycle — Wizs + Was = —(Qc + Qr) 
à 


Rendement Weycle 


=e 
Qc 


Qc 





en = | 


thermique du cycle 
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8.7.4.3 Cycle avec fractionnement de la compression 
Le principe est d’effectuer un refroidissement intermédiaire entre plusieurs 


compressionsadiabatiques de taux de compression moindre 


Cr 











W 


Cycle 











Figure 8.15 : Cycle de Joule avec fractionnement de la compression 


Remarque : 


Le Cycle de Joule comportant les trois améliorations à la fois est donné par la figure ci-dessus : 


Régéeneérateur 


© 
7 b 


D 
@ 

— Turbine I DD lr S 

Y re Woycle 





Retroidisseur 
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Figure 8.16 : Cycle de Joule avec régénération et fractionnement 
de la détente et de la compression 


8.8 Le Cycle de Stirling et d’Ericsson 

Les cycles de Stirling et d’Ericsson correspondent à des machines quasi mythiques que 
nous ne savons pas vraiment construire de façon opérationnellement utile. Ces cycles ont 
permis d'introduire certains concepts importants (exp: régénération) qui ont permis 


d’améliorer le rendement de vraies machines. 


=" Le cycle de Stirling se compose de 2 transformations isothermes (1-2), (3-4) et de 2 


transformations isochores (2-3), (4-1) : 





Figure 8.17 Cycle de Stirling 
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=" Le cycle d’Ericsson se compose de 2 transformations isothermes (1-2), (3-4) et de 2 


transformations isobares (2-3), (4-1) : 





Figure 8.18 Cycle d’Ericsson 


La régénération est dite parfaite si la totalité de la quantité de chaleur Q23 est récupérée 
totalement le long de la détente (4-1). 


Le rendement thermique du cycle de Stirling et d’Ericsson sont tels que : 


Nth Stirling — Ith Ericsson — 1+ QF/Qe = 1 — T;/Tc 
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Chapitre 9 : 


Changement de phase des corps purs 


9.1 Définition de phase et changement de phase 


= On appelle phase d’un système toute partie homogène (ou état physique) en équilibre. Les 
paramètres de l’équilibre seront la pression et la température. Il y a quatre phases de la 


matière : a)-Solide  b)-Liquide c)-Gaz d)-Plasma 





























Figure 9.1 Trois états de la matière 


Remarque : 


Le plasma est un état de la matière plus dilué que le gaz. Il est caractérisé par un désordre plus 
grand que celui du gaz, puisque même les éléctrons et les ions ne sont plus liés mais se 


déplacent librement les uns par rapport aux autres. 


= Le passage d’une phase a une autre (Fig. 9.2) s'appelle changement de phase (ou 
changement d’état) qui s’acompagne, pour un corps pur, d’un échange de chaleur sans 


variation de pression et de température : Sa deg V 4 Oe Vv 


fea 
A à | Figure 9.2 : Les différents 


| tion _ | changements de phase de 


la matière 
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Expérience : 


Cette expérience à pour but de montrer le changement de phase d’un gaz réel, chose qui 
ne peut être observé pour un gaz parfait. Pour cela on comprimé isothermiquement (à 
Tambiante) UN gaz pur (CO; par exemple) dans un cylindre. Si le gaz est parfait, il va suivre la 
loi PV=nRT=Cst et on devrait observer une hyperbole dans le diagramme (P,V) (voir 
Fig.9.3). En fait l’expérience nous montre qu’il n’en est pas vraiment ainsi : en comprimant 
isothermiquement le CO», on observe des paliers de pression. Le lèr palier correspond à la 
transision vapeur-liquide (liquéfaction), pour laquelle la compression ne modifie plus la 
pression du gaz (Ps=Cst) mais son état. Une fois le gaz devenu entièrement liquide, la 
compression permet a nouveau d’augmenter la pression du liquide (V ~ Cst, liquide est 
presque incompressible) jusqu’à l’obtention du 2°" palier de pression où l’on assiste à une 
autre transition solide-liquide (solidification). On observe le même phénomène si on détent 
isothermiquement le CO; initialement a l’état solide. Dans ce cas on passera d’abord par le 
ème 


1% Palier (vaporisation). 


pamier (fusion) puis le 


leme Paber ( Solidification ) 








(a) Vapeur (b) Liquide-vapeur (c) liquide 


[| Vapeur 
[ ] liquide 





Figure 9.3 Changement de phase d’un gaz réel lors d’une compression 


ou détente isotherme 
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9.2 Diagrammes d’état des gaz réels 


9.2.1 Diagramme (P, 7, V) 


De même que pour le gaz parfait, on peut voir dans l’espace (P,V,T) la surface 
caractéristique d’un gaz réel donnée par la figure 9.4 : 


Cr - Point critique 
Tr : Point triple 





Figure 9.4 Surface caractéristique d’un gaz réel 


La projection de cette surface permet d’avoir plusieurs diagrammes avec lesquels on travaille 


le plus souvent (diagrammes (P,7) et (P, V)). 


9.2.2 Diagramme (P,7) 

Le diagramme d’équilibre (Fig. 9.5) est composé de 3 branches qui délimitent les trois 
phases de la matières : 

= Branche BA : courbe d’équilibre S-V 

= Branche BD : courbe d’équilibre S-L, la pente de cette branche est positive ou négative 


= Branche BC : courbe d’équilibre L-V 
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Le point triple est un point pour lequel les trois phases de la matière S-L-V existent en même 


temps. 


Pia Substance qui Substance qui 
augmente de diminue de volume 
arden D pendant 
x volume pendant | Pendant sa 
LA sa congélation | congélation 
vA 
` | C Point cnhiqu 
\ 3 
3 % = 
e $ LIQUIDE 

AS F: 

o L .d 

LT 5 7 = 

k j rs 

a | 4% 
SOLIDE % | 
k ; 
B Point triple 
VAPEUR 
se 
got 
n i 
T: 
P [bar] 
10° | 


10° une an € 
solidification ee 
Liquide 7 Jiguéfaction 
10! Solide 
Fa Vapeur 
sublimation 
- s 


10 | ye 





A n 
-80 





-60 -40 -20 0 20 £40 T°c 
Figure 9.5 Diagramme de phase d’un corps pur dans le plan (TP) 
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9.2.3 Diagramme (P, V) 


Sur le diagramme (P,V) on observe les courbes de saturation solide-liquide, solide-vapeur 


et liquide-vapeur. La courbe de saturation qui nous intéresse le plus dans ce chapitre est la 


coube de saturation Liquide-vapeur, qui est sous forme d’une cloche dans la figure 9.6. Plus 


de détalis, concernant cette courbe, seront donnés ultérieurement. 







Cr Point critique 





LIQUIDE 


SOLIDE 


VAPEUR 


SOLIDE + LIQUIDE 










LIQUIDE + VAPEUR 


ligne triple 


SOLIDE + VAPEUR 


Figure 9.6 Courbes de saturation 


9.3 Courbe de saturation Liquide-Vapeur 


Soit un corps pur à l’état initial (i), si on l’échauffe à pression constante ce corps passe par 
les étapes suivantes (voir Fig. 9.7): 


i: liquide comprimé ou liquide sous refroidi 


a: le début de la vaporisation (changement de phase) qui se manifeste par l’apparition de la 


première 


bulle de vapeur. Le corps est à l’état de liquide saturé 


i-a : échauffement isobare du liquide 


b: melange 


liquide-vapeur saturé qu’on appelle aussi vapeur humide. Le liquide est en 


équilibre avec la vapeur. 


liquide. 


la fin de la vaporisation qui se manifeste par la disparition de la dernière goutte de 
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a-c : changement de phase liquide-vapeur qui se fait à pression et température constante (T7) 
f: vapeur surchauffée ou vapeur sèche 


c-f: échauffement de la vapeur à pression constante . 





Liquide comprimé 
Vapeur humide 


Liquide saturé 
Vapeur saturée 


vapeur surchau fée 


= 
k> 
D 

kr 


| Point critique 
i ot ¥ C 








I L a! k b i = ‘= f 
Liquide = mms à 
sn h 2%. Vapeur ~, 
comprimé | = AF N NE pn 
I = . Mélange + -surchauffée Li one — 
-4 Liquide- % apeur Le . = d; 
TR Saturé DEA 
, ‘Fo 
l 7 bs "ie, Sy 
. 


Liquide Vapeur 





Figure 9.7 : La courbe de saturation Liquide-Vapeur 


La courbe de saturation est composée de la courbe d’ébullition et de la courbe de rosée. La 
courbe d’ébullition est l’ensemble de points à l’état de liquide saturé et la courbe de rosée est 
l’ensemble de points à l’état de vapeur saturée. A l’intérieur de la courbe de saturation, on a 
un mélange liquide-vapeur en équilibre, à sa gauche on a du liquide comprimé et à sa droite 
(en dessous de l’isotherme critique Teritique) de la vapeur surchauffee. En dessus de l’isotherme 


critique, on a du fluide supercritique. 
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9.4 Pression et température de saturation 
La pression et la température de saturation (ou de vapeur saturante) sont la pression et la 
température auxquelles se produit le changement de phase. Elles sont, respectivement, notées 


par Psat et Tsat . Ces deux variables sont liées entre elles par la relation de Rankine: 








In Prat = Aq at 
N Psat = A1 — 
Tsat 
A, et B; sont des constantes propres à chaque gaz 
Psat 
| D LL LLL C „Poin cnitique 
| 
| 
| 
VAPEUR | 
Pr _ | 
* Point tre 
| | 
| | 
| | 
| | Figure 9.8 : La courbe de saturation 
T T Tan Liquide-Vapeur (P,T) 


Remarque : 


= Le changement de de phase L <—V est une transformation réversible. Pour pouvoir la 
réaliser, 1l faut traverser la courbe de saturation. Autrement dit, T< T< Ter et Pix P < Po 
( voir Fig. 9.8) 


= La pression de saturation est une fonction croissante de la température de saturation. 


9.5 Chaleur latente de vaporisation 

C’est la quantité de chaleur nécessaire pour faire passer une unité de la matière de l’état 
liquide à l’état vapeur (et vice versa). Cette chaleur ne provoque pas une augmentation de la 
température mais sert à changer la phase. La chaleur latente est noté par Ly , elle est fonction 


de la température selon la relation de Clapeyron : 


dPsat 


Ly = Tsat(v"’ — V') dT 





avec 
Ly : la chaleur latente de vaporisation, [kJ/kg] 


v" : le volume massique de la vapeur saturée, [m°/kg] (Fig. a) 
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v' : le volume massique du liquide saturé, [m°/kg] (Fig. a) 


AP;q,/aT : la pente de la courbe de vaporisation au point considéré (Fig. b) 





P | Pat 
Per = 
| 
| 
| 
| 
P | 
B dP, (aT | 
Ru 
| — Ï 
Tr Top sm 
y! y" y [ m'kg ] 
Figure a Figure b 


Dans le diagramme entropique (Fig. 9.9), la surface en dessous de la transformation qui 


caractérise le changement de phase (a-c) est équivalente à la chaleur latente : 
Ly = Tiat(S” — S") >0 (vaporisation, a-c) 
Ly = Tsat(s’ — S”) <0 (condensation, b-a) 


Remarque : 


La chaleur latente de vaporisation est nulle au point critique: Ly(T,) = 0, et elle est 
maximale au point triple : Ly(T.) = Ly max: 








a a 
Liquide ,~—> 
4 | _ vapeur 


Pas 
| Figure 9.9: La chaleur latente de 
s' s" s[kl’keK]  vaporisation Lv sur le diagramme entropique 
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9.6 Le titre massique de la vapeur humide (ou saturée) 


Le titre de la vapeur est une proprièté du mélange liquide-vapeur (c'est-à-dire de la vapeur 
humide) qui désigne le taux de vapeur dans ce mélange. On le défini comme étant le rapport 


de la masse de vapeur dans le mélange sur la masse totale du mélange : 


My masse vapeur 
x = — = C, 
My +m, masse totale 
= Cas de liquide saturé: my =0 > x=0 (courbe d’ébullition) 


= Cas de vapeur saturée: m, = 0 > x=1 (courbe de rosée) 


On peut aussi calculer le titre de la vapeur par la formule suivante : 


ab v-v 
X= 
ac 


yp" —y' 
v': le volume massique du liquide saturé = Vz / mz, 


v” : le volume massique de la vapeur saturée = Vy / my 


v :le volume massique du mélange = Vr/ mr 


vapeur 





v [m kg] 











Démonstration : 
L V 
Vr = Vy + V; => MTV = mv’ + myv” => V = í + p” > 
Mr Mr 
1 11 V = V 
v=(1-x)v+xv > x==~ 
V — v 


De façon analogue on peut montrer que le titre de la vapeur est égal à : 


h—h' s= u—u' 
x =—— = = 


h" —h’ En sl —s! 7 uw’ —y! 
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où: 
h',s',u’: enthalpie, entropie et énergie interne massique du liquide saturé, respectivement 


h”,s”,u” : enthalpie, entropie et énergie interne massique de la vapeur saturée, respectivement 


h, s,u : enthalpie, entropie et énergie interne massique du mélange liquide-vapeur, respectivement 


9.7 Les courbes isotitres : 
Une isotitre est le lieu des points de titre égale. Toutes les isotitres se rencontrent toutes au 


point critique Cr. L’allure de ces courbes sur le diagramme entropique (T,s) et le diagramme 


de Mollier (h,s) est donnée par la Figure 9.10. 


h 
pP = constant 


UE consan 


fe COSI ory 
ffs CONS tan 


m 
S 
|- 
Tl 
[= 
m E 
© 
il 





Diagramme de Mollier (41,5) 


Diagramme entropique (T5) 


Figure 9.10 Allure des isotitres dans les diagrammes (T,s) et (h,s) 
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Chapitre 10: 


Cycles de centrales thermiques à vapeur 


10.1 Introduction 

Le présent chapitre concerne les systèmes générateurs de puissance à partir de la vapeur 
d'eau (fluide moteur) produite par une source d'énergie telle que: la combustion d'un 
combustible fossile, l’ énergie solaire, l'énergie nucléaire.....etc. (ces systèmes sont la base de 


fonctionnement des centrales thermiques électriques à vapeur d'eau). 


10.2 Principe de fonctionnement général d'une centrale thermique à vapeur 

Les composantes fondamentales d'une centrale à vapeur sont montrées schématiquement 
par la figure ci-dessous. Le système global de la centrale peut être décomposée en quatre 
sous-systèmes À, B, C et D, que nous allons décrire succinctement. Notons que notre étude 
sera concentrée essentiellement sur le sous-système A, où s'effectue la conversion de 
l’énergie thermique (chaleur) en énergie mécanique (travail). 
Sous-système B: Fournit la chaleur nécessaire pour vaporiser l'eau qui passe par la chaudière. 
Dans les centrales à combustible fossile, ceci est réalisé par transfert thermique entre les gaz 
chauds produits par la combustion et l'eau qui traverse les tubes de la chaudière. Dans les 
centrales nucléaires, l'origine de la chaleur est une réaction nucléaire contrôlée ayant lieu dans 
un réacteur nucléaire isolé. 
Sous-systeme A: La vapeur produite dans la chaudière traverse une turbine, où elle se détend 
et sa pression baisse en faisant tourner l'arbre de la turbine (production de travail) branché à 
un générateur électrique (sous-système D). La vapeur sortant de la turbine traverse le 
condensateur, où elle se condense par contact avec des tubes acheminant l'eau de 
refroidissement. Elle traverse ensuite une pompe pour être de nouveau envoyée dans la 
chaudière où elle sera vaporisée. 


Sous-système C : Comporte le circuit de l'eau de refroidissement. 
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| 
| 
| 
| 
: |- Cheminée 
| 
| 
| 
| 
| 
| 
| 





| 
chaudière |, 
lévaporateur}e 





J 


Pompe eau de refroidissement 





Figure 10.1 : Schéma général d’une centrale thermique à vapeur 


10.4. Les cycles théoriques d'une centrale thermique 

Les connaissances de base nécessaires pour faire une analyse thermodynamique des cycles 
de centrales thermique incluent le principe de conservation de la matière, le principe de 
conservation de l’énergie pour les systèmes ouverts (1% principe de la thermodynamique) et le 
deuxième principe de la thermodynamique. Ces principes s'appliquent aux différents 
composants de la centrale: turbine, pompe, et échangeurs de chaleur (chaudière et 
condenseur). Nous allons également supposer que les composants de la centrale fonctionnent 
en régime permanent (débit du fluide moteur constant) et la variation de l'énergie cinétique 


(AE) et potentielle (AF,) négligeables. 


10.4.1. Cycle idéal de Rankine 

Ce cycle est représenté sur la Figure 10.2, il est composé des transformations réversibles 
suivantes : 
1-25 : détente isentropique dans la turbine ; 
2s—3: échange de chaleur (rejet) à pression constante dans le condenseur (ou condensation) ; 
3—4s: compression isentropique dans la pompe ; 
4s—1 : échange de chaleur (apport) a pression constante dans la chaudière avec 4s-a 


préchauffage isobare de liquide et a-1 vaporisation de la vapeur. 
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Figure 10.2 : Cycle idéal de Rankine sur le diagramme (T-s) 


La vapeur saturée à l’état (1) entre dans la turbine où elle se détend isentropiquement et 


produit du travail (W, ou puissance mécanique WT) : la pression et la température de la 


vapeur diminuent durant ce processus, à la fin duquel la vapeur se retrouve à l’état humide 
(2s). La vapeur humide traverse ensuite le condenseur où elle rejette de la chaleur 
(Qpr, ou puissance thermique Q f) dans une source froide (lac, rivière, atmospheére,...etc.) et 
se transforme en liquide saturé (3). L’eau a l’état (3) traverse la pompe où elle subit une 
compression isentropique jusqu’à la pression de service de la chaudière (notons que sur la 
Figure 2.10 la distance entre le point 3 et 4s est exagérée et ce pour des raisons de clarté du 
schéma). L’eau en sortant de la pompe, à l’état de liquide comprimé (4s), entre dans la 
chaudière où elle reçoit de la chaleur (Qç, ou puissance thermique Qc) à partir d’une source 
chaude et se transforme en vapeur saturée (1) qui se dirige vers la turbine afin de compléter le 


cycle. 


Remarque : Le cycle idéal de Rankine inclue aussi la possibilité de surchauffer la vapeur 
saturée ala sortie de la chaudiere (c.a.d avant son entrée dans la turbine) comme le montre 
le cycle 1’-2’s-3-4s-1’ de la Figure 2.10. La surchauffe (1-1’) se fait dans un appareil appelé 
surchauffeur, et l’ensemble chaudière - surchauffeur porte le nom de générateur de vapeur. 
L’avantage de cette opération est de réduire la teneur de la vapeur humide en liquide (x2,, > 
X25) qui risque d’endommager (par érrosion) les aubes de la turbine. De ce fait, cet appareil 


devient moins efficace et necessite plus de maintenance. 


Résumons toutes les étapes du cycle idéale de Rankine (avec et sans surchauffe) dans le 
tableau suivant puis effectuons son bilan: 
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a) Les différentes étapes du cycle 
1-25 Détente isentropique dans la turbine (production de travail) : Q 2=0 
Wis = Wi, = m Ahy2, = Ms — hy) < 0 , [kW] 
|'-2’s 
( cycle avec ) 
surchauffe 


Wi = are =m As — Mm(hz,s — hy,) , [kW] 


45-3 Condensation de la vapeur (chaleur rejetée à P et T Cste) : W>,3 = 0 
Qs =m Ahzs3 = Mhz — hs) < 0 , [kW] 


re 
( 


Qr = m Ahz,s3 = M(h3 — hr) , [kW] 

surchauffe 

Compression isentropique dans la pompe (Q 34=0) 

Wi E Was = m Ahz4, = Mas — hz) > 0, [kW] 

WE =m i v dP ~ mvz3(P,, — P3) (fluide incompressible v = cst) 


Préchauffage (4-a) et vaporisation (a-1) isobares dans la chaudière : War = 0 


Qc = mMAhgs, = MR — has) > O , [KW] 
4s-1' 
( cycle | Qc =m Ahgsi, = m(hy = has) > [kW] 
surchauffe 





m désigne le débit massique de l’eau (fluide moteur) dans la centrale, [kg/s] 


Remarque : Sur le diagramme (T-s) (Fig.10.2), la quantité de chaleur reçue par la source 
chaude est équivalente à la surface (c-4s-a-1-b) en dessous de la transformation 4s-1 (ou 
encore la surface (c-3-a-1’-b') en dessous de la transformation 4s-1' pour le cycle avec 
surchauffe). La quantité de chaleur cédée à la source froide est équivalente a la surface (b-2s- 
3-c) en dessous de la transformation 2s-3 (ou encore la surface (b'-2’s-3-c) en dessous de la 
transformation 2’s-3 pour le cycle avec surchauffe). La somme algébrique de ces quantités de 
chaleur est égale au travail du cycle équivalent a la surface (3-4s-a-1-2s-3) (ou la surface (3- 


4s-a-1’-2’s-3) pour le cycle avec surchauffe). 


b) Le bilan du cycle 





. . : oo LOT OP : ; 
Puissance utile Wautite — W cycle — Wis T Wis J —(Qc + Qr) 
es” 
oT 
Whet 


Rendement W — hə. 
thermique du cycle — |- — = ko 
S 
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( cycle avec ) 


surchauffe 





WT, est la puissance nette de la turbine, égale à la somme algébrique de la puissance brute 


(WT) délivrée par la turbine et la puissance nécessaire à la pompe (WP). 
is P P pomp 


c) Effets de la pression de vaporisation et de condensation sur le cycle de Rankine 
Considérons la transformation isobare 1-2 évoluant entre l’entropie Sı et S2 (Fig. 10.3). La 
chaleur Q;2 mise en jeu dans ce cas correspond bien à la surface en dessous de la 
transformation en question. La transformation isotherme a-b évoluant entre S; et S2 tel que 
Qap = Q12 (surface S,-a-b-Sp = surface S;-1-2-S.) est équivalente à 1-2. En partant de ce 
principe le cycle a-b-c-d-a (qui est un cycle de Carnot) est équivalent au cycle de Rankine 
]-2s-3-4s-5-1. La relation qui s’établie dans ce cas entre les rendements thermodynamique des 


deux cycles est la suivante : 


Tp 
(Mtn)rankine = (Mtn)carnot équiv — 1 — T 
a 


Ta et T1} sont des températures moyennes de la source chaude et froide du cycle de Carnot 


équivalent 









Carnot 
équivalent 





51-54 S3e Sh , 


Figure 10.3 : Cycle de Carnot équivalent au cycle de Rankine 


D’après la relation ci-dessus, le rendement thermique du cycle (tp) Rankine augmente si la 
température moyenne de la source chaude (T,) augmente ou que celle de la source froide (Tp) 
diminue. De plus, comme T, et Tp varient respectivement de la même façon que la pression 
de vaporisation (P1=P;= P4s) et de condensation (P2,= P3), nous en déduisons que le 
rendement (7:n) rankine augmente si la pression de vaporisation (dans la chaudière) augmente 


(Fig. 10.4a) ou que la pression de condensation (dans le condenseur) diminue (Fig. 10.4b). 
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Augmentation 
de la pression T 


de vaporisation Pression de 
vaporisation 
fixée 


Diminution de 
} atm | 








la pression de 


Pression de condensation 


condensation 


- 00° C 
fixée 100” C 


f 
B 
2b 


Ambient 
temperature 





h) 5 


(a) (b) 


Figure 10.4 : Amélioration du rendement thermique du cycle idéal de Rankine. 
(a) par augmentation de la pression de vaporisation. (b) par diminution de la pression de 
condensation. 


Remarque : D’après la Figure 10.4a, lorsqu'on augmente la pression de vaporisation pour 
améliorer le rendement thermique du cycle on fait croitre en même temps l’humidité de la 
vapeur à la fin de la détente! (X24 < X2) = présence plus importante de gouttelettes d’eau 
dans la turbine = usure mécanique par érosion des ailettes = réduction du rendement global 
de la turbine. Pour éviter cette situation nous appliquons une resurchauffe intermédiaire dont 
le principe est présenté dans le paragraphe suivant. De la même façon, en diminuant la 
pression dans le condenseur (Figure 10.4b) on augmente l humidité de la vapeur à la fin de la 


détente (X2p < x3). 


10.4.2. Cycle à resurchauffe 

Dans ce cas la détente isentropique se fait dans deux étages de la turbine (transformations 
isentropiques 1-2s et 3-4s, Fig. 10.5) entre lesquelles la vapeur subit une resurchauffe 
(transformation isobare 2s-3). La resurchauffe accroit le taux de la vapeur dans la turbine 


(X45 > Xas) et augmente de 4 à 5% le rendement thermique du cycle’. 


a) Les différentes étapes du cycle 


Détente isentropique dans le 1“ étage de la turbine : Q ;,=0 
W, = Wiz2s =m Ah = m(hzs — hı) < 0 , [kW] 


Resurchauffe dans le resurchauffeur (apport de chaleur à P Cste) : W,3 = 0 
Qcz = M Ahzs3 = Mhz — hs) > 0 , [kW] 





* Cette situation a lieu même dans le cas du cycle idéal de Rankine avec surchauffe. 

? On peut multiplier le nombre de resurchauffe pour augmenter davantage le rendement du cycle. Cependant, 
au-delà de deux resurchauffes cette augmentation et assez faible pour justifier la complexité et le coût de 
Pinstallation. 
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D 


Détente isentropique dans le étage de la turbine : Q34=0 


W,.? = Wea, = m Alga, = Mas — hz) < 0 , [KW] 


Condensation de la vapeur (chaleur rejetée à P et T Cste) : Wis = 0 
Qf =m Ahass = m(he — has) < 0 : [kW] 

Compression isentropique dans la pompe (Q56=0) 

Wis = Wees = th Ahses = Res — hs) > 0 , [kW] 

WE =m Je v dP ~ mv.(P., — Ps.) (fluide incompressible v = cst) 


Production de la vapeur dans le générateur de vapeur (apport de chaleur à 
P=cst): We, = 0 


Qc, = m Ales, = M(Ngs — h1) > 0 , [kW] 





b) Le bilan du cycle 


—— 
. | Wret 


Rendement Weycie — Qf 14+ hs — has 
thermique du cycle Qc 








Oc (Aes — h) + (hz — hs) 














eau de f 
refroidissement 


Figure 10.5 : Cycle de turbine à vapeur avec une seule resurchauffe 


10.4.3. Cycle à régénération 
Dans ce paragraphe nous allons étudier l’évolution d’un cycle à régénération assez simple, 


car en pratique il existence d’autres versions plus complexes. 
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La particularité du cycle de régénération consiste à détendre la vapeur dans le premier étage 
de la turbine (1-25) à la sortie duquel une fraction (y) du flux total est soutirée et envoyée dans 
un mélangeur’. La fraction restante de vapeur (1-y) poursuit sa détente isentropique dans le 
deuxième étage de la turbine (2s-3s) puis traverse le condenseur pour se transformer en 
liquide saturé (3s-4) qui sera comprimé par la pompel (4-5s) jusqu’à la pression de service du 
mélangeur (P2;). On mixe dans le mélangeur (y) quantité de vapeur à l’état (2s) avec (1-y) 
quantité de liquide à l’état (5s) pour obtenir du liquide saturé (6) qui sera comprimé par la 
pompe? (6-7s) jusqu’à la pression de service du générateur de vapeur (Pı). Pour finir, le fluide 


est chauffé de l’état (7s) à l’état (1) dans le générateur de vapeur. 





Figure 10.6 : Cycle à régénération à un seul mélangeur 


Remarque : sur la figure 10.6, le point (a) représente l’état du fluide à la sortie de la pompe 


pour un cycle sans soutirage (dans ce cas là Qe = h4 — ha ) 


Comment déterminer les fractions de vapeur soutirées à chaque étage de la turbine ? 
La conservation du débit massique à travers un volume de contrôle renfermant les deux 
étages de la turbine (Fig. 10.7) se traduit par la relation suivante : 
M3 + M = M 


En divisant par m, on obtient: 


3 7 XX s X 4 aN A . 2 . 2 
C’est un échangeur de chaleur à contact direct où deux fluides à la même pression et de température différentes 
sont mélangés pour former un fluide à une température intermédiaire. 
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M m 
—+—=1 

My M 
Si on note par y (= m,/m,) la fraction de vapeur soutirée au 1“ étage de la turbine (point 2), 


la fraction du flux total soutirée au 2°" étage de la turbine (point 3) est telle que : 


Figure 10.7 : Bilan de matiere pour un volume 
de controle entourant la turbine 





De la même façon, la conservation de l’énergie” pour un volume de contrôle entourant le 
mélangeur supposé parfait: pas d'échange de chaleur avec l'extérieur (Fig. 10.8), permet 
d’écrire : 
0 = mhs + Maihs — mihe 
En divisant par m; on obtient : 
0 = yh + (1 — y)hs — he 
ou encore : 
he — hs 
Y = hah; 
Donc sachant les états du fluide au point (2), (5) et (6), on peut déterminer la fraction de 


vapeur (y) soutirée au 1“ étage de la turbine et celle soutirée au second étage (1-y). 
p y £ g y 


Figure 10.8 : Bilan d’énergie pour un volume 
de contrôle entourant le mélangeur 





* Il s’agit du 1° principe de la thermodynamique écrit pour un système ouvert à deux entrées et une sortie, en 
régime permanent et avec une variation négligeable de l’énergie cinétique et potentielle (voir Chapitre 2). 
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a) Les différentes étapes du cycle 


Détente isentropique dans le 1°” étage de la turbine : Q 2=0 
W, = Wigs = m Ahizs = Thy (hzs — hy) <0 , [kW] 


LS 


Détente isentropique dans le 2°™ étage de la turbine : Q23=0 


W? = W535 = m3 Ah2535 = M3 (hss — hs) < 0 , [kW] 
Condensation de la vapeur (chaleur rejetée à P et T Cste) : W34 = 0 
Qs = M3 Ahgs4 = M3Ch4 — has) < 0 , [kW] 


Compression isentropique dans la pompel (Q45=0) 


y 1 E Wiss = M3 Auss = M3 (hss — hy) > 0, [kW] 


Wr. = m; Í H v dP ~ m3 V,(Ps, — P,) (fluide incompressible v = cst) 


Mélange de (y) vapeur avec (1-y) liquide dans le mélangeur 


0 = Mhs + Mhs =a Mhe 


Compression isentropique dans la pompe2 (Q57=0) 
W 57 = Wéss = m4 Ahezs = (hg — hzs) > 0, [kW] 


We =m, J 4 v dP = m, ve(Pzs — P) (fluide incompressible v = cst) 
Production de la vapeur dans le générateur de vapeur (apport de chaleur à 
P=cst) : Wo, = 0 


Qc = Mm, Ahzs1 = Mu (hys — hy) > 0 , [kW] 





b) Le bilan du cycle 


Puissance utile Wie = Voie =WI14+WI2 + we + W,.? = —(Qc + QF) 
_—_——— SS LL 
T 


We 





Rendement ; mz (h4 — 3s) 
thermique du cycle : : mu (h-s — hy) 
S 


10.5. Influence des irréversibilités et des pertes sur les cycles de vapeur 
Les irréversibilités et les pertes sont associées à chaque élément de la centrale à vapeur. 
Certains de leurs effets ont une influence plus marquée sur le rendement thermique du cycle 


que d'autres, particulièrement dans la turbine et la pompe. 
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10.5.1. Cas de la turbine 

La présence des irréversibilités lors de la détente adiabatique (1-2) dans la turbine 
s'accompagne d'une augmentation d'entropie (Fig. 10.9). La puissance mécanique WT 
obtenue dans ce cas (réel) est inférieure à la puissance mécanique WT obtenue dans le cas 
idéal (1-2s détente adiabatique réversible). Le rapport de ces puissances est dit rendement 
isentropique (ou mécanique) de la turbine : 


_ puissance réel W" h;—-h 
hos p hy 


Nis 


puissance idéal Wī 


P] 








Détente réel 
(irréversible) 


/ Détente idéal 


/ (réversible) p, 
j 4 - 


Cai 





(a) (b) 


Figure 10.9 : Représentation de la détente adiabatique réversible (idéal) 
et irréversible (réel) sur le diagramme T-s (a) et h-s (b) 


Comme les irréversibilités diminuent la puissance de la centrale (Wycie), par réduction de la 


puissance brute de la turbine (WT < W% ), elles réduisent le rendement thermique du 
cycle qui s’exprime par : 
\Weycte| = [WT + Wi 2 Nis Wi + Wi 


He Ge Oc Oc 


10.5.2. Cas de la pompe 

La aussi les irréversibilités causent une augmentation d'entropie (Fig. 10.10). Mais la 
puissance mécanique W? obtenue dans le cas réel (3-4: compression irréversible) est 
supérieure à la puissance mécanique W£ obtenue dans le cas idéal (3-4s : compression 
réversible). Le rapport de ces puissances est dit rendement isentropique (ou mécanique) de la 
pompe : 


p puissance idéal Wy, has — hz 





flis puissance réel WP h,—hg 
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h P, 4 
Ri 
ha. Compression 
4s | - Se ee 
veel (irrév) 
Compression 


idéal (réver) 


P; 
2 





Figure 10.10 : Représentation de la compression adiabatique réversible (idéal) 
et irréversible (réel) sur le diagramme T-s (a) et h-s (b) 


Les irréversibilités dans la pompe diminuent la puissance de la centrale (Weycte); par 
augmentation de la puissance de la pompe (W? > WÈ ), ce qui réduit le rendement thermique 
du cycle : 

ae [Weyetel _ [ws + WP ul WP] _ | Wis + Wis /nis” | 

Qc Qc Qc 
Remarque : si on considère les irréversibilités dans la turbine et la pompe a la fois, le 
rendement thermique du cycle est tel que : 
|Weyete| _ [WT + WP| _ | nis Wis + Wis / nis” | 


eG: Oc Oc 
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Une machine frigorifique est un appareil destiné a créer du froid. Pour cela, il doit enlever 
(en un cycle) une quantité de chaleur Q; à la source froide (exemple: air intérieur du 


réfrigérateur) et donner une quantité de chaleur Qc à la source chaude (exemple : air ambiant 


de la cuisine), tout en recevant un travail W, du milieu extérieur, (Fig.11.1). 


Qc <0 Qf =0 





SYSTEME 


FLUIDE 
THERMIQUE 





TRAVAIL 
W >0 


Figure 11.1 Schéma de fonctionnement d’une machine frigorifique 


Les machines frigorifiques sont des machines réceptrices (reçoivent du travail) qui se divisent 


en deux grandes classes, selon l’ organigramme présenté ci-dessous : 


Machines frigorifiques 


à compression à absorption 
de vapeur de gaz 
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Contrairement aux machines frigorifiques a absorption et ceux a compression de gaz, les 
machines frigorifiques a compression de vapeur sont les appareils les plus répondus. Toutes 
ces machines fonctionnent selon des cycles qui utilisent des fluides appelés fluides 
frigorifiques ou réfrigérants. Dans le cas des machines frigorifiques a compression de vapeur, 
le réfrigérant les plus utilisés depuis de nombreuses années est le dichlorodifluorométhane 
(CFC de nom commercial R12). Depuis 2001, il est interdit d’utiliser les CFC qui ont été 
remplacés par le R134a. Ce réfrigérant est également utilisé avec les machines a compression 


de gaz. Pour les machines à absorption le réfrigérant est une solution à base d’ammoniac. 


L’efficacité thermique d’un cycle de machine frigorifique est le rapport entre l’énergie 
utile et l’ énergie coûteuse qu'il a fallu fournir. Autrement dit, le rapport de ce qu’on veut (Qp 


sur ce qu'il faut dépenser pour l’avoir (W). 


énergie utile _ Q; 


énergie côuteuse Weycle 





Eth >ou<1 





Tous les cycles des machines frigorifiques sont parcourus dans le sens contraire des 


aiguilles d’une montre (ou sens trigonométrique). 


Le cycle de Carnot est théorique il est le premier cycle à être imaginé pour faire 
fonctionner une machine frigorifique. Il est composé de deux transformations isothermes 


(2-3) et (4-1), et de deux transformations isentropiques (1-2) et (3-4). 


Compression isentropique du réfrigérant de l’état vapeur humide (1) à l’état vapeur 


saturée (2s) (O2, =0) : 
Wis = W 495 = Ahy25 = hs — hy > 0 
Condensation : le réfrigérant cède (a P et T constantes) la chaleur Qc à la source 


chaude et passe de l’état de vapeur saturée (2s) à l’état de liquide saturée (3) : 


Qc = Tc(s3 — S25) < 0 
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Détente isentropique du réfrigérant : 


WE — W 34, = Ah34, = Ras = hz <0 


Vaporisation : le réfrigérant enlève (à P et T constantes) la chaleur Q; à la source 
froide (production de froid) et se vaporise partiellement de l’état (4) à l’état (1). Le 


titre de la vapeur x; > x4: 


Qr = T; (sy — S45) = Tr (S25 —s3) >0 





ci al a 
B o T 






Efficacité 


thermique du cycle 


Source chaude 
a Te 


Source froide 


‘1 





Figure 11.2 Cycle frigorifique de Carnot 
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"= L'efficacité du cycle frigorifique de Carnot représente la valeur maximale théorique de 


tout cycle frigorifique qui opère entre deux températures Tc et T; 


= Les deux cas pour lesquels l'efficacité thermique d’un cycle frigorifique est > 1 ou < 1 


sont illustrés par la figure suivante : 


2s 





4s. U d 
NAN , 











= 


= Le cycle de Carnot évolue entièrement a l’intérieur de la courbe de saturation, ce qui 


présente les inconvénients suivants : 


> le titre de la vapeur (x;) avant son entrée dans le compresseur est incertain. Il y a 
donc risque que la chaleur de compression soit insuffisante pour vaporiser tout le 
liquide aspiré. Il en résulte une diminution de la capacité d’aspiration du 


compresseur. 


> on risque également les coups de liquide lorsque le volume du liquide restant en fin 


de compression est supérieur au volume mort (V2). 


> le film d’huile qui lubrifie les surfaces internes du cylindre ne devient plus efficace 


car il est lavé par le liquide restant. 
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Le cycle idéal (Fig. 11.3) diffère du cycle de Carnot par la présence d’une vanne de détente 
à la place de la turbine et par le fait que ce cycle n’évolue pas entièrement dans la zone de la 


vapeur humide. 


Miheu extérieur chaud 


P o- 













4e 









À 


' Vapeur sat 





(A) @) 


Figure 11.3 Cycle idéal d’une machine frigorifique représenté sur 
(A) diagramme (T-s) — (B) diagramme (P-h) 
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Compression isentropique du réfrigérant de l’état vapeur saturée (1) à l’état vapeur 
surchauffée (2s) (Q)2; =0) : 
W‘, = Wiz = Ahi; = hos = hı >0 


Refroidissement isobare puis condensation : le réfrigérant cède la chaleur Qc à la 
source chaude et passe de l’état vapeur surchauffée (2) à l’état liquide saturée (3) : 

Qc = Ahzs3 = h; — hz, < 0 

Détente isenthalpique (de Joule Thomson) du réfrigérant à travers la vanne de détente 
(Q34 = 0 et W34 = 0) : 

Ah3, = 0 > h, = h, 

Vaporisation : le réfrigérant enlève (à P et T constantes) la chaleur Q; à la source 
froide (production de froid) et se vaporise de l’état de vapeur humide (4) à l’état de 
vapeur saturée (1). 


Qr = Ahy, = h -h4 > 0 





Qy hı — h4 


Efficacité 


. Eth =, 
thermique du cycle W cycle 





= Toutes les transformations du cycle idéal sont réversibles sauf la transformation 


isenthalpique (3-4) qui est irréversible. 


=" Dans le cycle idéal d’une machine frigorifique (Fig.11.3), le fluide frigorifique subit une 
détente isenthalpique (3-4) et non pas isentropique (3-4s) comme dans le cas du cycle de 
Carnot. Cela entraine une perte de l'effet frigorifique (surface (4s-4-S4-S4-4s), Fig.11.3A 
ou segment (4-4s) Fig.11.3B). Pour éviter cette perte, on doit refroidir le fluide frigorifique 
à la sortie du condenseur de façon à récupérer la quantité de froid perdue. Cette question 


fera l’objet du prochain paragraphe. 
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Afin de récupérer la quantité de froid perdue à cause de la détente isenthalpique, on 
refroidit le liquide à la sortie du condenseur par la vapeur qui sort de l’évaporateur, à travers 


un échangeur de chaleur auxiliaire. De ce fait, le cycle actuel d’une machine frigorifique est 


donné par la figure 11.4. 


Milieu extérieur chaud 


Fo 








Compresseur 
Echangeur W, 
de chaleur ve 
auxiliaire 


Vanne de at - 
détente A Ur 





Espace refregire froid 











$a'= 33 
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a) différentes étapes du cycle 


Compression isentropique du réfrigérant à l’état de vapeur surchauffée (Q72; =0) : 


Wi, = Wins = Ah — hs = hi, > 0 


Refroidissement isobare puis condensation : le réfrigérant cède la chaleur Qc a la 
source chaude et passe de l’état vapeur surchauffée (2s) à l’état liquide saturée (3) : 
Qc = Ahzs3 = hz — hg, < 0 
Refroidissement isobare du réfrigérant (liquide) dans l’échangeur de chaleur 
auxiliaire : 

Qsr = Ahg3, = hx — hz <0 

Détente isenthalpique (de Joule Thomson) du réfrigérant a travers la vanne de détente 
(Q34 = 0 et W3-4=0): 

Ah3,4 =0 > hy = hy 

Vaporisation : le refrigerant enlève (a P et T constantes) la chaleur Qs à la source 
froide (production de froid) et se vaporise de l’état de vapeur humide (4') à l’état de 
vapeur saturée (1). 

Qr = Ahzan = hy — hy > 0 

Echauffement isobare du réfrigérant (vapeur surchauffée) dans l’échangeur 
auxiliaire : 

Qsr = Alyy, = hy — hy > 0 (= h; —h,’) 





b) Le bilan du cycle 
Energie coûteuse | Woycie = Wig = —(Q¢ + Qc) 


Efficacité Qe hı = hy 


R = — e O 
thermique du cycle |" Weycie hy — Ay + hzs—h3 





11.3.3 Cycle cascade de réfrigération 

Certaines applications industrielles nécessitent de basses températures qui ne peuvent être 
atteintes par un seul cycle de réfrigération. Pour cela on utilise plusieurs cycles de 
réfrigération qui opèrent en série. Un tel cycle est appelé cycle cascade de réfrigération. La 
figure 11.5 représente un cycle cascade de réfrigération à deux étages (I et IT). L’ échangeur 
de chaleur intermédiaire joue le rôle d’évaporateur pour le cycle I et le rôle de condenseur 


pour le cycle II. 
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la capacité de 
| réfrigération 
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Milieu extérieur chaud 


he. 


© 





Vanne de 
\ détente 





(étage IT) 





Milieu réfrégiré froid 


Figure 11.5: Cycle cascade de réfrigération à deux étages 


Remarque : 


Le débit du fluide frigorifique dans chaque 


même (Mm; € My,[Kg/s]). Si on suppose que 


étage du cycle n’est pas obligatoirement le 


le transfert thermique dans l’échangeur de 


chaleur du centre, est parfait, on peut écrire la relation suivante : 


My (hs = hg) = —My] (hz — has) 
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Chapitre 11 


Le bilan du cycle 


Weycle = (WS), F (WS), = m,(h25 — hy) + my (hes — hs) 


Puissance coûteuse 
= -(Q, + Qc) >0 


Efficacité Q, Ch — h4) 
thermique du cycle |“? ~ W cycte ho — hy) + Ty (hes — hs) 





11.3.4 Cycle de réfrigération à compression multi-étagé 

Lorsque les fluides frigorifiques utilisés dans les différents étages de la cascade (Fig.11.5) 
sont de même nature, on remplace l’échangeur de chaleur intermédiaire (entre le cycle I et II) 
par une chambre de flash (voir Fig.11.6) qui permet de séparer les deux phases liquide-vapeur 
du réfrigérant a la sortie de la vanne de détente (6) : une quantité de liquide saturée (1-x6) est 
envoyée vers la deuxième vanne de détente (7) et une quantité (xs) de vapeur saturée est 
envoyée vers une chambre de mélange (3). Par la suite ces deux quantités (c'est-à-dire (1-x6) 


et x6) vont se mélanger dans cette chambre. 


Le bilan du cycle 


Op = (1 —6)(hy — hg) > 0, [kJ/kg] 


Energie coûteuse | Weycte = (Wis), + (Wis), = (1 — %6) Chas — h1) + Chas — ho) > 0 


Efficacité Q f (1 — x6) (Ch; — hg) 
Se — aa 
thermique du cycle i Weyce (1 — x6)(hzs — h1) + (has — ho) 





Le titre de la vapeur x6 est calculé a partir de la relation suivante : 


_h—h 
hg —hy 





X6 
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Chapitre 11 


L’enthalpie (hg) est déterminée a partir du bilan d’ énergie de la chambre de mélange : 
énergie à l'entrée = énergie à la sortie 


Xéh3 + (1 — xs)hzs = (1)ho 


Mibeu extérieur chaud 


# (rc 









45 
travail de il i 
i hy # ! ! H, 
. compression “s, Se i 
T | FE + Vanne de ns 
4505 + détente 


Compresseur 
$to haute pression 






chambre de 
flash 






chambre de 






la capacité de | 
réfnigérahon Vanne de 


détente 


(l-X 6) 








ee 





basse pression 





Mibeu réfregire froid 


Figure 11.6 Cycle de réfrigération a compression multi étagé 


11.4 Cycle de réfrigération à compression de gaz 


L’ensemble des cycles frigorifiques considérés ci-dessus nécessite un changement de phase 
du réfrigérant. Par contre pour le cycle de réfrigération à compression de gaz, le réfrigérant 
reste toujours à l’état gazeux. Ce cycle permet d’atteindre de très basses températures 
(-150°C) nécessaires à la liquéfaction des gaz et d’autres applications particulières. Le cycle 
de réfrigération à compression de gaz est le cycle inverse de Brayton (parcouru dans le sens 
contraire des aiguilles d’une montre) (voir Fig. 11.7). Il est composé de deux transformations 
isobares (2-3) et (4-1), et de deux transformations isentropiques (1-2s), (3-4s). Les 


transformations (1-2) et (3-4) sont simplement adiabatiques (non réversibles). 
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Niheu chaud 






| Compresseur 


Echangeur 
de chaleur 


Cr 


Miheu Froid réfrigere 


Figure 11.7 : Schéma d’une machine frigorifique a compression de gaz 





Figure 11.8 : Cycle de réfrigération à compression de gaz 


(a) cycle réversible - (b) cycle irréversible 
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a) Bilan du cycle réversible 


Energie rejetée Qc = h; — hv, < 0, [kJ/kg] 


Energie coûteuse Weycle = Wis + Wis = (Rs a hy) + (has B h3) >0 


Efficacité Q f hy — has 
CO — E 
thermique du cycle VO Weycte (has — hy) + (has — h3) 





Bilan du cycle irréversible 
Energie perdue Qc = h; — h, < 0,[k]/kg] 


Energie utile Qs = h — h, > 0, [kJ/kg] 


Energie coûteuse Weycle =W° +W?” = (h2 — hy) + (hy — hg) 


Efficacité Q hı — h4 


i = P S E E ___ 
thermique du cycle |"  Weycie (hz — hi) + (hg — ha) 
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Annexe 


Travail de transvasement 
Soit un fluide de masse m présent à la frontière (1) d’un système ouvert (voir Fig.a). Pour 
introduire cette masse dans le système, il suffit de la déplacer d’une distance Ax,. Le travail 


Wir, nécessaire pour cela (lorsque la force F; exercée est constante) est égal a: 
V4 
Wer = F,. Ax, = Fe — Pj: Vi 


P; et V, désignent, respectivement, la pression et le volume de la masse m à la frontière (1) 
De façon analogue (et en tenant compte de la convention des signes en thermodynamique), le 


travail W,,., nécessaire pour faire sortir la masse m du système ouvert, est égal à : 


V2 
Wir, = —F,. Ax» = ae = —P;.V, 


P; et V, sont, respectivement, la pression et le volume de la masse m à la frontière (2) 


Le travail de transvasement total W}, est alors exprimé comme suit : 


W; = Wa + We, = P,.V, — P2. V» 





n. Système ouvert a 
frontiere frontière 


(1) (2) 


Figure a : Travail de transvasement 


Identités analytiques 


Soit une fonction à trois variables indépendantes f (x, y,z) = 0. La dérivée totale de fest : 


af af Of 
a ca En 0 


Si l’on se place a z constant (dz = 0), il reste 





Of Of 
x ox + a = 0 
soit : 

(a5) 
dy _ \ax 
= & 

dy 


: dy p . Oy $ gax 
Comme z est constant le quotient => représente aussi ~~ Où l’on considère que x est une 


fonction de y et z. On peut donc écrire : 





ap 
03.5- j 
On aurait de même 
ap 


(5). | = 


En faisant le produit (*) x (**), on obtient 


= 
[Ss 
ST 


(=) = = (xxx) 


De façon analogue à (*), on peut écrire qu’à x constant puis à y constant : 


(3z) 


D). (22) 








aE 


T (2) (xxx) 


Le produit membre à membre des égalités (** ) et (**** ) nous permet d’écrire : 


G) G Go), =-1 
Oy! \ðz/y \Ok)y.- 


Ce résultat peut en effet être généralisé à une fonction a n variables f (x1, X2, ...,X,) = 0, tel 


que : 
OX: OX) OXn 
a pag ees ieee — =f n 
(=) i Gr) ‘ toy ae 


X2, X3» Xn—-1 


